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ABSTRACT 


Heat transfer performance of horizontal roped tubes was 
studied using an evaporator/condenser test rig. The condenser 
section contained a single horizontal tube and was operated at 
both vacuum and atmospheric conditions. A vacuum pump was 
integral to the system to ensure the presence of 
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experimentally obtaining the overall heat transfer 


coefficient, U and then utilizing a modified Wilson Plot 
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Results obtained for steam condensation indicate that the 
MHT Korodense titanium tube provides negligible enhancement 
compared to a titanium tube. The addition of wire-wrapping 
provided minimal increase in the outside heat transfer 
Seerficient, about 8% at a fractional coverage of 
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Korodense and Wire-Wrapped Smooth Tubes at 
Atmospheric Pressure 

Frictional Temperature Rise for the MHT 
Korodense Titanium Tube with and without a 


HEATEX Insert 
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NOMENCLATURE 
as defined in equation (4.30) 


effective anside surftacewencea mime 


ef fective outside surtacemeane amma 

as defined in equation (4.26) 

as defined in equation (2.22) 

specific heat at a constant pressure, J/kg-K 
as defined in equation (2.14) 

leading coefficient for the inside heat transfer 
correlation used 

as defined in equation (2.8) 

as defined in equation (2.23) 

external groove diameter, m 

tube inside diameter, m 

tube outside diameter, m 

wire diameter, m 

groove depth, m 


as defined in equation (2.32) 


as defined in equation (2.33) 
as defined in equation (2.35) 
as defined in equation (2.30) 
as defined in equation (2.4) 
gravitational constant, 9.81 m/s 


specific enthalpy of vaporization, J/kg 


inside heat transfer coefficient, W/m -K 


Stiecnidle eat sttugamiamets “Cocwma1 cleiieee W/m - K 

outside heat transfer coefficient for enhanced tubes, 
W/m? -K 

outside heat transfer coefficient for smooth tubes, 
W/m--K 

thermal conductivity of the coolant film, W/m-K 
thermal conductivity of the condensate film, W/m-K 
thermal conductivity of the tube material, W/m-K 

as defined in equation (2.26) 

as defined in equation (2.27) 


active length of tube exposed to steam, m 
external groove depth, m 

length of inlet portion of tube, m 

length of outlet portion of tube, m 

log mean temperature difference, K 
Reynolds number exponent 

mass flow rate of the coolant, kg/s 
as defined in equation (4.25) 
Prandtl number exponent 

number of tubes in a vertical row 
Nusselt number 

COuruGaeloen pitch, m 

wire-wrap pitch, m 


Saturation pressure, Pa 


Prandtl number 


X1 


Pr 


PW 


Prandtl number at the wall temperature 
as defined in equation (2.710) 
heat transfer rate, W 


heat. fice mn 


tube specific constant as defined for equation (2.30) 


radius of curvature of corrugation ridge, m 
fouling resistance, K/W 
inside thermal resistance, K/W 
outside thermal resistance, K/W 
total thermal resistance, K/W 
wall thermal resistance, K/W 

Reynolds number 

Reynolds number for the condensate film 


two phase Reynolds number 


temperature difference across the condensate film, K 


temperature difference across the film, K 
mean coolant bulk temperature, K 

vapor saturation temperature, K 

outside wall temperature, K 

cooling water inlet temperature, K 

cooling water outlet temperature, K 
overall heat transfer coefficient, W/m*-K 
vapor velocity, m/s 

coolant velocity, m/s 

Reynolds number exponent in equation (4.15) 


as defined in equation (4.24) 


eae 


as defined in equation (4.23) 
as defined in equation (4.19 or 4.20) 


as defined in equation (2.9) 


Greek Symbols 


Oe 


Pe 


Py 


leading coefficient for the outside heat transfer 
correlation used 

helix angle measured from the tube axis 

as defined in equation (2.37) 

as defined in equation (2.25) 

enhancement ratio based on AT for wire-wrapped smooth 
tubes 

enhancement ratio based on AT for wire-wrapped 
corrugated tubes 

as defined in equation (2.29) 

mass flow rate of condensate per unit length of 
horizontal condenser tube 

dynamic viscosity of the coolant, kg/m-s 

dynamic viscosity of the condensate film, kg/m-s 
dynamic viscosity of condensate at the wall, kg/m-s 
surface efficiency 

as defined in equation (2.14) 

tube specific constant as defined for equation (2.10) 
as defined in equation (4.15, 4.16 or 4.17) 

density of the condensate film, kg/m° 


density of the vapor, kg/m? 


Bese 


surface tension of the condensate at film temperature, 
Lbe/ Be 
severity factor 


tube specific constant as defined for equation (279) 


Dea 


ACKNOWLEDGEMENTS 

The author would like to express his sincere appreciation 
Memeroressor Paul J. Marto and Professor Stephen B. Memory for 
their advice, guidance, and continual support towards the 
completion of this thesis. 

Many thanks is also given to Mr. Tom McCord and the 
workers in the Mechanical Engineering Department Machine Shop 
for their timely fabrication of all the required materials 
necessary to complete this endeavor. 

The author wishes to extend his greatest appreciation and 
thanks to his wife, Laura. Without her extreme patience, 
love, and encouragement this project would not have been 


possible. 





i. 2v-RODUELEON 


A. BACKGROUND 

Reductions in size, weight and cost have become 
increasingly important design considerations aba the 
development of every component destined for use aboard United 
States Navy vessels. Significant reductions in weight have 
already been achieved for main propulsion machinery with the 
meeene Of the gas turbine propulsion plant. Weight reduction 
is also possible through redesign of auxiliary systems and 
other large components such as main steam condensers. The use 
of heat transfer enhancement techniques for condenser tubing 
has been one avenue receiving considerable attention in the 
pursuit of this goal. Smaller more efficient condensers would 
result in reductions in capital cost, space requirements, 
weight and operating cost. 

The four primary factors that affect the performance of 
condensers are the thermal resistances on the coolant side, 
across the tube wall, the vapor side, and also that due to 
fouling. Each can be reduced to minimize its effect, which in 
turn enhances the overall heat transfer rate of the system. 
In experimental studies, the fouling resistance can be 


Megiected if the tube is thoroughly clean prior to testing. 


Enhancement techniques can be categorized into two major 
groups, active and pasSive. Guido [Ref. 1] compiled a fairly 
comprehensive listing of both groups of enhancement techniques 
for the tube side heat transfer rate. Passive techniques such 
aS internal helical ribbing and displaced promoters (the use 
of wire mesh or HEATEX inserts) will be discussed in detail in 
a later section. It should be noted, Nowever, that both of 
these techniques have a corresponding increase in pumping 
power requirements. Wall thermal resistance can be reduced 
passively by decreasing the wall thickness or changing the 
tube material. Similar to enhancement techniques for the tube 
Side, the shell side resistance can also be reduced by both 
active and passive means. The passive enhancement techniques 
of wire-wrapping and tube corrugation will also be covered in 
a later section. 

A variety of corrugated or roped tubes are available for 
use in condenser applications. Figure 1 shows the profile of 
one type of corrugated tube which has a groove radius smaller 
than the ridge radius. Wolverine Tube Company manufactures a 
low pressure drop (LPD) and a maximum heat transfer (MHT) tube 
of this type (trade name: Korodense), available for commercial 
use. The MHT tubes have a larger groove depth compared to the 
LPD tubes [Ref. 2]. Other types of corrugated tubes areyaee 
available which have a larger groove radius than the ridge 
radius. Marto et al. [Ref. 3] present test results on various 


other corrugation geometries. 








Figure 1. Profile of a corrugated tube 


B. CONDENSATION APPLICATIONS FOR ENHANCED TUBES 

In 1971 Withers and Young [{Ref. 4] performed a theoretical 
analysis on a power plant main condenser. Their findings 
indicated a substantial saving in almost every condenser 
aspect if the smooth tubes were replaced by Korodense roped 
tubes. Similar results were obtained by Webb et al. [Ref. 5]. 
They indicated a possible 30 to 40 percent increase in the 
Overall heat transfer coefficient when analyzing a tube-for- 
tube replacement of smooth to Korodense tubes. This benefit 
was calculated for two separate power plants of different heat 


duty. An ongoing comprehensive study of 14 Tennessee Valley 


AUCHOrI Cy sete power plant condensers, retubed with 
corrugated tubes, has yielded significant economic and 
performance benefits [Refs. 6 and 7]. Of the 14 condensers 
retubed by TVA, their Gallatin Unit 1 has provided long-term 
results with over 12 years of successful operation since 
retubing. In 1983 Boyd et al. [Ref. 8] studied the 
performance of TVA’s Gallatin Unit 1 plant after only three 
years of operation. Their findings yielded up to age: 
increase in the overall heat transfer coefficient when clean 
tubes were analyzed. However, data also indicated that the 
corrugated tubes fouled quicker than plain tubes. Rabas et 
al. [Ref. 9] studied 12 of the retubed TVA condensers and 
found fouling rates for the enhanced tubes ranging from about 
the same as to almost twice that experienced in plain tubes. 
Although the fouling rates were higher, Rabas et al. [Ref. 9] 
reported the fouling resistance values for the LPD Korodense 
tubes never exceeded the minimum allowable TEMA standard. 
Union Carbide [Ref. 10] performed analysis on roped 
condenser tubes with a Thin Film Condensing Promoter 
(T.F.C.P.), which is essentially a wire-wrapped roped tube. 
Their findings indicated an increase in the heat transfer 
coefficient by a factor of between 2 to 3 times that found 


with a smooth tube under condensate loaded conditions. 


C. NAVAL POSTGRADUATE SCHOOL CONDENSATION RESEARCH 

Derantumeemas a mumber Of properties which make it a 
desirable alternative to the conventional copper-nickel tubing 
currently in use today on board Naval vessels. Bes high 
strength to weight ratio enables the tubes to be designed with 
a much thinner wall thickness, while still maintaining the 
durability of the overall system. Titanium also has a high 
resistance to erosion and sea water corrosion. This leads to 
reduced shutdown time for tube cleaning and replacement. One 
of the primary disadvantages of titanium tubing is its poor 
Boaermal conductivity. Titanium has approximately 50% the 
Miermmal Conductivity of copper-nickel [Ref. 11]. This 
disadvantage can be partially offset by a decreased wall 
thickness and also through the use of enhanced surfaces. The 
Naval Postgraduate School (NPS), in conjunction with the Naval 
Surface Warfare Center, Carderock Division, Annapolis 
Detachment, has been testing the enhancement characteristics 
of various types of titanium tubing. 

In order to examine the heat transfer characteristics of 
various types of enhanced tubes, a suitable test platform must 
be utilized. While the ideal test platform would be a full 
Scale model of the actual condenser used in the field, the 
cost associated with such an apparatus would be prohibitive. 
The condensation test platform, constructed by Krohn [Ref. 12] 
at NPS, has been in use since 1982. Although several 


modifications have been made by various researchers since that 


time, the overall apparatus and the operation remain basically 
unchanged. Van Petten [Ref. 13] lists a comprehensive summary 
of the work completed by eight NPS researchers from 1982 
Bhivewgh. 1 ore Most of the work concentrated Cm 
performance of various types of externally finned euvess 
Mitrou [Ref. 14], in addition to studying finned tubes, also 
performed experiments on wire-wrapped smooth copper tubes. 
More recent work conducted by Coumes [Ref. 15], Guttendorf 
[Ref. 16] and Swensen [Ref. 17] dealt with further studies of 
finned copper tubes. O’Keefe [Ref. 18] studied the 
performance of various titanium tubes. His tests included 
smooth tubes, wire-wrapped smooth tubes, LPD Korodense tubes, 
and wire-wrapped LPD Korodense tubes. Other researchers 
including Kanakis, Brower and Mazzone [Refs. 19-21] have 
studied tube bundle configurations for smooth and corrugated 


tubes, both with and without wire-wrapping. 


D. OBJECTIVES 
The main objectives of this study are as follows: 


1. To verify the repeatability oof the calculated 
coefficients for both coolant side and vapor side heat 
transfer correlations obtained by previous researchers on 
wire-wrapped smooth and corrugated titanium tubes. 


2. To manufacture and collect condensation data on smooth 
and MHT Korodense titanium tubes with varying fractional 
coverage due to wire-wrapping. 


3. To determine if any relationship exists for optimal 
fractional coverage or pitch to diameter ratio between 
wire-wrapped smooth tubes and wire-wrapped MHT and LPD 
Korodense tubes. 





Il. LITERATURE SURVEY 


A. INTRODUCTION 

Condensation on horizontal tubes has been an intensely 
Studied heat transfer mechanism throughout the years. New 
developments pal condenser tubing require cont inued 
experimental and theoretical research in this field. 

There are two fundamental modes of condensation on 
horizontally mounted condenser tubes, filmwise and dropwise. 
Dropwise condensation can produce heat transfer rates as much 
as an order of magnitude larger than those associated with 
filmwise condensation. Various surface coatings have been 
developed in order to promote dropwise condensation and take 
advantage of the higher heat transfer rates. However, with 
time these coatings lose their effectiveness and filmwise 
condensation eventually returns. For this reason, condenser 
design calculations are normally based on the more 
conservative filmwise mode of condensation [Ref. 22]. 

Internally, the primary factors affecting thermal 
resistance are the intensity of turbulence and the thermal 
conductivity of the coolant medium. The degree of turbulence 
can be changed by altering the internal tube wall geometry or 


through the use of displacement promoters. 


B. VAPOR SIDE CONSIDERATIONS 

Condensation on horizontal tubes consists of a complex two 
phase process. A suitable theoretical model must be utilized 
in order to accurately predict these phenomena. 

1. Horizontal Smooth Tubes 

Nusselt [Ref. 23] developed the foundation for the 

study of filmwise condensation on horizontal tubes in 1916. 
His work was formulated for a quiescent vapor condensing ona 
Single horizontal tube. Due to the increase in the thickness 
of condensate as gravity forces it around the sides of the 
tube, the local heat transfer coefficient decreases around the 
tube circumference. Nusselt developed a theory for the 
average heat transfer coefficient around the tube, accounting 
for the lower resistance at the top of the tube where the film 
was the thinnest and the higher resistance at the bottom of 
the tube where the film was the thickest. As a result, he 


developed the following equation: 


1/4 
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where the fluid properties are evaluated at the film 


temperature, given by: 
7, =r ee (2.2) 
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Nusselt’s theory has been extensively studied since 


1916, and with the imposed assumptions, has been found to be 





generally valid [Refs. 24, 25]. It has also been found to be 
meeurate £or cases which do not conform to Nusselt’s original 
assumptions, such as variable wall temperature [Ref. 26]. One 
of the major problems encountered in applying Nusselt’s theory 
in the design of condensers arises from his assumption of a 
quiescent vapor. While in theory, and in some limited 
practical applications, the assumption of a stationary vapor 
can be justified, most steam condensers operate under 
conditions where the vapor is travelling at some constant 
velocity. 

With the vapor in motion, sufficient shear forces can 
be developed which significantly thin the condensate film and 
increase the outside heat transfer coefficient. Barly 
theoretical work by Shekriladze and Gomelauri [Ref. 27], which 
accounted for vapor shear, resulted in the following equation 


to approximate the mean Nusselt number. 
= 0.64(1 +(1 + 1.69F) 1/2) 1/2 (2.3) 


where: 
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and the two phase Reynolds number is given by: 


pe. Ugis 
H 


R€,4 = (2.5) 


Lee and Rose (Ref. 28] compared vapor shear models with 
experimental results and found that the Shekriladze-Gomelauri 
results were more conservative than more rigorous developments 


performed by other researchers due to their simplified 


approximation for the interfacial shear stress. Fujii Geeaie 
[Ref. 29], in a more recent study, performed an extensive 
development on the outside of a horizontal tube. From 
experimental data, they also developed an empirical 


formulation for the condensation of steam on a horizontal tube 


which includes vapor velocity effects: 
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where F and Reo, are defined in equations (2.4)and (2.5). For 
Situations where the surface shear forces dominate, equation 
(2.6) tends to more accurately predict the vapor side heat 
transfer coefficient for steam. 

As discussed previously, the primary modes of external 
tube enhancement covered in this study are wire-wrapping and 
the addition of helical grooves. Both methods have associated 
advantages and disadvantages. 

2. Wire-Wrapped Smooth Tubes 

The benefits of wire-wrapping smooth tubes have been 

attributed to two primary effects: the thinning Of Ep 


condensate film between the wires due to surface tension 


forces and the better drainage of the condensate along the 


10 








wire, (Fujii et al. [Ref. 30] and Marto et al. [Ref. 31]). 
Thomas [Ref. 32] conducted tests on vertical tubes with wires 
attached vertically along the tubes and demonstrated this 
latter enhancement effect. Horizontal smooth tubes with 
coiled wire were later studied by Thomas et al. [Ref. 33] and 
their findings indicated a two fold increase in the heat 
transfer coefficient over that predicted by Nusselt theory for 
a smooth tube. Marto et al. [Ref. 31] determined an optimum 
wire pitch to diameter ratio between 5 and 7 for steam 
condensation. They achieved an enhancement between 1.5 and 
1.8 for the heat transfer coefficient over plain tubes. 

For low surface tension fluids (R-11 and ethanol), 
Fujii et al. [Ref. 30] obtained results 2 to 3 times those 
predicted by Nusselt theory, with an optimum wire pitch to 
diameter ratio of 2. O'Keefe [Ref. 18] compared the 
enhancement ratio versus fractional coverage (the percentage 
of the tube’s surface covered by wire) on a titanium smooth 
tube. Together with data obtained by Mitrou [Ref. 14] and 
Sethumadhavan and Rao [Ref. 34], he determined that the 
optimum fractional coverage was between 10 and 30%. This 
fairly broad range of values was attributed to the possible 
effects of the different tube material thermal conductivities 


used in the studies. 


Atak 


3. Corrugated Tubes 
In 1971, Withers and Young [Reft:. 35] perieeieeaw. 
design study of two condenser bundles with heat duty of 25 
million Btu/hr, one with bare (smooth) tubes and the other 
with corrugated (roped) tubes. Coolant flow pressure drop was 
one of the matching parameters for the comparison. This 


requirement made it necessary to use a larger corrugated tube 


(1.00 inch outside diameter (OD)), than the smooth tube (0.75 
inch OD). They found an enhancement of 88% for the outside 
bundle heat transfer coefficient. Catchpole and Drew [Ref. 


36] tested single corrugated tubes and determined that the 
steam side heat transfer coefficient increased with increasing 
groove density and passed through a maximum with increasing 
groove depth. Their conclusion was that the smaller the helix 
angle (measured from the tube axis), the worse the condensate 
drainage, resulting in a reduction of the outside heat 
transfer coefficient. Similar results were reported by Marto 
et al. [Ref. 3] and Mehta and Rao [Ref. 37]. Eleven different 
corrugated tubes of three primary types were tested by Marto 
et al. [Refw 3]. Within each corrugation type, the amgess 
value for the outside heat transfer coefficient was obtained 
at small pitch to diameter ratios (P/D). As P/D increased, 
the enhancement ratio decreased. In tests conducted by Mehta 
and Rao [Ref. 37], they found the condensing coefficient 
enhancement ratio decreased from 1.26 to 1.1 when eae 


corrugation groove pitch was increased from 6.35 mm to 12.5 mm 


a2 





with a constant groove depth of 0.63 mm. They also tested the 
merece, Of holding the groove pitch to groove depth ratio 
constant. A maximum enhancement ratio of 1.38 was achieved 
with the smallest groove depth and pitch. Zimparov et al. 
[Ref. 38] confirmed the same trend, using 27 mm outside 
diameter corrugated tubes, with the enhancement ratio 
Piemedasing from 0.99 to 1.22 for a corrugation groove pitch 
decreasing from 16.9 mm to 6.5 mm. 

Cunningham and Bahjernejad [Ref. 39] examined the 
effects of condensate inundation on corrugated tube bundles. 
They tested two different multistart tubes, one with two 
helical starts and the other with six Starts. Their results 
indicated severe limitations in performance due to inundation 
effects for tube bundles more than a few rows’ deep. 
Additional inundation tests were conducted on multistart tubes 
by Cunningham and Boudinar [Ref. 40]. They reported similar 
results as those of Cunningham and Bahjernejad [Ref. 39], and 
additionally found that the more severe the indentation 
(groove depth) on the roped tube, the greater was the 
reduction in vapor side heat transfer for the tube bundle. 

Rabas [Ref. 41] studied an extensive data bank of 
results published by numerous researchers. He found that the 
enhancement in the outside heat transfer coefficient was 


negligible for tubes with grooved geometries that are suitable 


aS 


for ball cleaning! of the inside "Surface. He also Comecmsed 
the collected data with the results predicted by three 


different correlations. The Wolverine correlation ({kKefgeZs 
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Ce = Za (2.8) 
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@® Values for w and @ are listed for specifi vom. 
geometries by Wolverine in Reference 2. 


@® The variable C, in Equation (2.7) must be applied using 
English units. 


and the correlation developed by Mehta and Rao [Ref. 37]: 
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tball cleaning is a method of removing deposits from the 
inside of condenser tubes by periodically recirculating sponge 
rubber balls through the tubes with the coolant. 
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moe 15x10 * < § < 200x10°*, were both found to overpredict the 
outside heat transfer coefficient. Values similar to those 
obtained from the experimental data were predicted by the 


Yorkshire correlation [Ref. 42] given by: 
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where: 


(2.14) 











wage Varies with tube profile and can be found in Reference 
42. However, Rabas [Ref. 41] cautioned its use due to the 
results being close to unity. In conclusion, he recommended 
using smooth tube correlations and neglecting the minimal 
enhancement associated with corrugated tubes for design 
calculations. 

Dreitser et al. [Ref. 43] determined a correlation for 
the enhancement ratio of corrugated tubes based on three 
primary geometric ratios: ijm——=hne Latae of the radius of 
curvature of the corrugation ridge to the outside diameter, 


Ry D 2) the ratio of groove pitch to the outside diameter, 


or 
P/D,; and 3) the ratio of the external groove diameter to the 


piecide diameter, ds/D.. 
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They found that if any of the three ratios were increased 
while keeping the others constant, the enhancement ratio would 
decrease. 

Cunningham and Fidler [Ref. 44] made a comparative 
study of the effects of vapor shear on multistart corrugated 
tubes and smooth tubes. They found that the corrugaledieames 
maintained their enhancement advantage over the smooth tube 
when subjected to increased vapor velocity. 

4. Wire-Wrapped Corrugated Tubes 

Wire-wrapping of corrugated tubes has been shown by 
various researchers to Significantly reduce the effect of 
condensate inundation (Memory et al. [Ref. 45], Union Carbide 
[Ref. 10], and Marto and Wanniarachchi [Ref. 46]). 

Union Carbide [Ref. 10] tested steam condensation for 
a wire-wrapped roped tube in a bundle arrangement. They 
reported nearly identical performance for both the un@eages 
and the condensate loaded conditions tested. Condensing 
coefficients approximately 2.5 times that of plain tube 
results were achieved for the condensate loaded condition. 
Memory et al. [Ref. 45] obtained a 60 to 68% enhancement over 
smooth tube results when condensing R-113 on a wire-wrapped 


roped tube bundle arrangement containing four tubes in a 


i 








vertical column. Their tests indicated an optimal wire pitch 
Mmemdlameter ratio of 7.5. 

O'Keefe [Ref. 18] tested the enhancement effects of 
three different wire diameters (1.6 mm, 1.0 mm, and 0.5 mm) 
when used to wrap an LPD Korodense tube. Enhancements ranging 
between 20 and 30% above smooth tube results were obtained 
with the maximum obtained for the 0.5 mm wire with a wire 


pitch to wire diameter ratio of 14. 


C. COOLANT SIDE CONSIDERATIONS 
1. Smooth Tubes 
Numerous correlations have been developed for use in 
the analysis of the coolant side heat transfer correlation 
associated with the turbulent flow regime inside a pipe (Re > 
10,000). Many of the earlier correlations have been presented 


inne form: 


Nu = C,;Re™Pr” (2.16) 


This form has been used for several well-known 
correlations including those developed by Dittus and Boelter 


[Ref. 47]: 
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and Colburn [Ref. 48] 
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Nu = 0.023 Re®?:*®§pri/3 (2.18) 


A correction factor for equation (2.18) was developed 


by Sieder and Tate [Ref. 49]: 


Nu = 01027 Reo? Pri) Ge ee (2.19) 


This modification compensates for the variation in the coolant 
viscosity when large temperature differences exist between the 
bulk coolant and the inner tube wall temperature. All fluid 
properties for equations (2.17), (2.18), and (2,05) 


evaluated at the mean coolant bulk Geriperael as 


qT, = is (2.20) 


In more recent years, Sleicher and Rouse [Ref. 50] and 
Petukhov [Ref. 51] developed equations (2.21) and (2.24) 
respectively, which are applicable over a wider range of 


Prandtl numbers. 


Nu = 5 + 0.015 Re,° Pr,,® (2 222s) 
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where: 
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Merenz et al. [Ref. 52] compared their experimental results 


fOr cold water (i.e. large Prandtl number) with seven of the 
more common inside heat transfer correlations. They found 
Ehat the Petukhov-Popov and Sleicher-Rouse correlations agreed 
to within 5% of the experimental data. The other correlations 
typically underpredicted the data by between 5 and 15%. 
2. Corrugated Tubes 

The internal ridges of the tube wall produced when 
manufacturing corrugated tubes cause two primary effects. The 
irregularity of the wall surface results in an increased 
frictional effect and causes the flow to swirl, thereby 
requiring additional pumping power to achieve flow rates equal 
to similar diameter smooth tube flow. The increase in 
turbulence also has the effect of increasing the inside heat 
meameter Fate. In the study by Withers and Young [Ref. 35], 
the total inside enhancement ratio of the tube bundle tested 
was 56% above that obtained for a smooth tube bundle designed 
with an equal pressure drop. Catchpole and Drew [Ref. 36] 
found the inside heat transfer coefficient to increase with 


increasing groove depth and groove frequency. However, they 


eo 


stressed that the increase in performance was obtained at the 
expense of significantly greater friction factors and therefor 
increased pumping power and operating cost. In comparing 
their results with the Prandtl analogy [Ref. 53] they found 
that, due to the turbulence caused by the presence of grooves, 
the inside heat transfer coefficient followed the trend 
predicted by the Prandtl analogy. 

Mehta and Rao [Ref. 37] presented their findings in 
terms of the severity factor, the ratio of groove depth 
Squared to the product of groove pitch and inside diameter. 
This dimensionless group was chosen because it has been found 
to be a useful criterion for evaluation of heat transfer and 
frictional characteristics of plain and spirally Gmmaameed 
tubes. An improvement of 260% was obtained above smooth tube 
results with the higher inside heat transfer coefficient 
corresponding to higher severity factors. Of the eleven tubes 
tested by Marto et al. [Ref. 3], inside enhancement was 
reported for all the tube geometries examined, with the 
magnitude ranging from 1.6 to 5.1 times a smooth tube. 
Zimparov et al. [Ref. 38] tested nine tubes with differing 
corrugation geometries and obtained results similar to 
previous researchers. They achieved inside enhancement ratios 
between 1.77 and 2.73, with an increase in the friction factor 
Ole Zoe ine leels 

For corrugated tubes, Withers ([Ref. 54] modified 


existing correlations and developed the following: 
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where: 

y= - 2.5 In2£] + 3.75} (2.29) 
D; 

and: 
SS ee ae a, ee (2.30) 
8 246 wry, a7 / re). | 


The variables m and r were determined for each tube by a 
Surve-fitting procedure. He obtained agreement between 
experimental and predicted data to within 10%. Wolverine 
currently utilizes a form of this correlation in their 
Engineering Data Book II [Ref. 2]. 

In a comparison of the LPD and MHT Korodense tubes 
manufactured by Wolverine [Ref. 2], they report a 25% higher 
inside heat transfer coefficient for the MHT Korodense tube 
than for the LPD Korodense tube. The MHT Korodense tube also 
has twice the pressure drop of the LPD Korodense tube. 

After compiling an extensive collection of published 
data, Rabas et al. [Ref. 55] recommended a minor modification 
to the leading coefficient in the Yorkshire correlation [Ref. 
42) to better fit the experimental data. The Yorkshire 


correlation is similar in form to equation (2.17) with an 


Pak 


additional complex leading coefficient (E,) to compen=agemaan 


the corrugation tube geometry. 


Nu =(0 30275. ege aa ae (2.31) 

where: 
Bosal, O + 1. 82 SBe Bey Gey red (25323 
EB, = GOs &) (22353) 
C= tame whe) (2.34) 
In(£,) = - © a,(1n(e/D;) 17 (2.358 
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elliGh: 


a, = 5.0766, a, = 4.4224, a, =1.1240, 42, = 6.700 ie 


The modification by Rabas et al. (Ref. 55] changed the 
exponent in equation (2.32) sfrom©.7 cameetoce 
Panchal and France [Ref. 56] conducted a study of 


large pitch spirally indented tubes for possible application 





to pressure drop sensitive systems. Their data resulted in 
internal enhancements as high as 100% over smooth tube data 
and they reported agreement within a 20% error band for the 
results predicted by the correlation proposed by Mehta and Rao 


[Ref. 57], given by: 
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— = (22-36) 
Remi? 2° Regt / 2 lake ey 
where: 
Bie Oia Ay (2 alae 1.2). ** (er37) 
and y is defined as in equation (2.29). It should be noted 


that 6 in equation (2.37) is constant for any given tube 
configuration. 
3. Tubes With Inserts 

Numerous types of inserts are available for use in 
enhancing the inside heat transfer coefficient. Some of the 
types used by previous researchers at NPS include: twisted 
tape, wire-wrap, and HEATEX (a wire mesh insert supplied by 
CAL GAVIN). A number of correlations have been developed for 
determining the inside heat transfer coefficient when using 
these various inserts [Refs. 58-60]. These correlations for 
HEATEX, however, are for specific loop densities of the wire 
matrix and are not universally valid. HEATEX inserts not only 
Significantly increase the overall heat transfer rates, but 
also inhibit fouling. Gough and Rogers [Ref. 61] presented 
results collected from a condenser fitted with HEATEX inserts 
which was in operation for over four years without being shut 
down for cleaning. Their findings indicated the actual 
fouling resistance was over 96% less than the allowable limits 


after a service life of 26 times that of the original unit 


ZS 


with no insert. However, it should be noted that the use of 
an insert has a pronounced effect on increasing the pressure 
drop across the tube. Inserts are used in experimental 
applications to realize a larger temperature drop in the 
coolant flow, thus decreasing the uncertainty of the heat 


Erancter ¢calethawrtencr 
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III. EXPERIMENTAL APPARATUS 


A. SYSTEM OVERVIEW 

The apparatus used for this research was originally 
constructed by Krohn [Ref. 12]. The initial design has been 
modified by several previous researchers (in particular 
Swensen [Ref. 17]) and has been used successfully by Coumes 
(Ref. 15], Guttendorf [Ref. 16], Swensen [Ref. 17] and O’ Keefe 
(Ref. 18] to test condensation on single horizontal tubes of 
various configurations. Figure 2 contains a general schematic 
of the overall system. The boiler, a 0.30 meter diameter 
Pyrex glass cylinder, contains ten 4 kW, 440 Volt Watlow 
immersion heaters which are used to generate steam from 
distilled water. From the boiler, steam passes up through a 
cylindrical section of Pyrex glass with an inside diameter of 
0.15 meters and a length of 2.13 meters. Two 90 degree Pyrex 
Glass elbows redirect the steam back down a second similar 
cylindrical section of 1.52 meters in length. The steam then 
enters the stainless steel test section containing the 
horizontal condenser tube mounted as shown in Figure 3. The 
test section also contains a circular glass viewing port so 
the condensation process can be observed while conducting 
tests. Excess steam, not condensed on the sample tube, passes 


through the test section and into the final Pyrex glass 
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Schematic of the Single Tube Test Apparatus 


Figure 2. 
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Figure 3. Schematic of the Test Section 
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cylinder containing an auxiliary condenser. The auxiliary 
condenser section collects all the condensate, where it 
returns to the boiler via a gravity drain in the baseplate. 

Cooling water for the single horizontal tube in the test 
section was supplied by two centrifugal pumps connected in 
series and drawing suction from an adjacent cooling water 
sump. A throttling valve and calibrated flowmeter allowed for 
control of the cooling water flow rate. The auxiliary 
condenser cooling water waS supplied directly from the 
building water main. After setting the desired coolant flow 
rate through the single horizontal tube, the Sdatumatien 
temperature in the test section was controlled by adjusting 
the throttle valve on the auxiliary condenser coolanmyimae7 
line. 

A suction port at the base of the auxiliary condenser 
section allowed for the removal of non-condensible gases. 
Figure 4 shows a schematic of the associated system. Non- 
condensibles drawn from the auxiliary condenser passed through 
an internal condensing coil where, any associated steam that 
might be present was condensed and collected ina plexiglass 
vacuum chamber for later removal. The non-condensibles then 
passed through the vacuum pump and were expelled to the 


abMmospaente. 
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Schematic of the Purging System and Cooling Water 


Figure 4. 
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B. SYSTEM INSTRUMENTATION 

The 440 Volt immersion heater power supply was controlled 
by a silicon controlled rectifier mounted on the electrical 
Switchboard. Poole [Ref. 62] details the power calculation 
for the input into the data acquisition system. 

Three separate methods were used to monitor the pressure 


of the system: 


1. Setra model 204 pressure transducer. 
2. System saturation temperature converted into pressure. 


3. Heise solid front pressure gage (visual reference only). 


Temperature measurements were obtained using several 
devices. The vapor temperature of the system was measured by 
a Teflon coated and a metal sheathed type-T copper/constantan 
thermocouple located just upstream of the horizontal tube 
being tested. Condensate return and ambient temperatures were 
also measured with Teflon coated type-T copper/constantan 
thermocouples. The inlet and outlet coolant temperature of 
the condenser test tube were measured with three separate 


probes: 


1. A Teflon coated type-T copper/constantan thermocouple. 


2. A ten junction, Teflon coated type-T copper/constantan 
thermopile. 


3. An HP 2804A quartz crystal thermometer. 
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The three temperature probes at the outlet of the horizontal 
tube were located downstream of a mixing chamber in order to 
obtaln a representative average coolant discharge temperature. 
Figure 3 contains the locations of these temperature probes. 

An HP-3497A data acquisition system, linked to a HP-9826A 
computer, was used to collect the voltage readings from all 
the measurement devices. The data collection and reduction 
program "DRPML", listed in Appendix C, allowed for the 
processing and storage of all the data by the HP-9826A 
computer. 

Swensen [Ref. 17] provides calibration data for the 
pressure transducer, thermocouples, quartz thermometer and 


flowmeter. 


Se SOUBES TESTED 

Eight enhanced tubes were manufactured for this thesis. 
Tests were also conducted on tubes studied by previous 
researchers to ensure repeatability of results. Table I 
contains a listing of all the tubes tested and additional 
tubes used in the comparison of results. All of the tubes 
assembled for this study were made from a single titanium MHT 
Korodense tube manufactured by Wolverine. Seven of the tubes 
were then wrapped with wire having diameters ranging from 0.1 
to 1.6 mm to test enhancement effects. To attach the wire to 
the tubes, two different methods were used. The larger 


diameter (stiffer) wires (2 0.5 mm) were wrapped around a rod 


eal 


of a smaller diameter than the tube, this formed the wire into 
a tight coil. The wire was then fed onto the tube by twisting 
the coil and advancing it along the groove of the tube. The 
wire was then spot welded in place at each end of the active 
condensing surface. For the wires < 0.5 mm, one end was spot 
welded to the tube at the edge of the active surface. The 
wire was then tightly wrapped around the tube, by hand, in the 
groove and spot welded to the other end of the active surface. 
Both methods resulted in a constant wire pitch equal to the 
corrugation pitch of the tube. The wire was not physically 
attached to the active tube surface and did not act as a fin. 
The additional tubes included in Table I were manufactured for 


research conducted by O’Keefe [Ref. 18] and Mitrou [Ref. 14]. 


a 


iS TINGSOR TUBES TESTED 


Table I. 
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IV. EXPERIMENTAL PROCEDURES AND DATA ANALYSIS 


A. TUBE PREPARATION 
To ensure consistent results in data collection, a 
continuous film of condensate must be present over the entire 
tube surface. A chemical treatment procedure used by 
Guttendorf [Ref. 16] and several other researchers at NPS was 
used to produce filmwise condensation on the tube surface. 
Prior to installation into the test section, the tube was 
chemically treated as follows: 
1. The inside and outside of the tube were thoroughly 


cleaned with mild soap using a soft bristle brush. The 
Soap was rinsed off with distilled water, then the tube 


was rinsed with acetone. Distilled water was again 
poured over the tube, to ensure the water formed a 
continuous film on the tube surface. If breaks in the 


film were observed (an indication of impurities on the 
tube surface) the cleaning process was repeated. Oncea 
continuous film was present, the active surface was not 
handled, thus ilaby dl akg awk ial the possibility of 
contamination. 


2. The tube was then placed in a steam bath. 


3. A solution of 50% by weight sodium hydroxide and ethyl 
alcohol was mixed for treating the tube surface. The 
solution was kept warm to ensure a watery consistency. 


4. The solution was then applied over the tube surface, with 
a small brush, in 10 minute intervals for one hour. 
While applying the solution, the tube was rotated so the 
entire surface could be treated. If the tube was treated 
previously, the solution was only applied every 5 minutes 
over a 20 minute period. 


5. After the tube surface was treated, it was removed from 


the steam bath and rinsed with distilled water to remove 
any excess alcohol/sodium hydroxide solution. To verify 
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a continuous film will be produced in the presence of 
Steam, the tube was returned to the steam bath, to ensure 
the entire surface wets easily. 


6. Without touching the active surface, the tube was removed 
from the steam bath and immediately mounted in the test 
section of the boiler apparatus. If tests were to be 
conducted with an insert, the insert was installed prior 
to securing the coolant flow piping to the test section. 


The oxide layer formed by the chemical treatment promotes 
a continuous film on the tube surface. Due to the negligible 
thickness of the oxide film, its thermal resistance was 
considered insignificant compared to the overall thermal 


resistance of the tube. 


B. SYSTEM START-UP AND SHUTDOWN PROCEDURES 
After the tube has been installed in the test section, the 
system can be brought on-line using the following procedures: 


1. The boiler section must be filled with distilled water to 
a level approximately 4 to 6 inches above the top of the 
heater elements. The boiler can be filled by attaching 
a hose from the distilled water tank to the piping at the 
boiler drain/fill valve. Opening the vent valve on the 
auxiliary condenser allows for gravity filling of the 
boiler. To drain the boiler, disconnect the distilled 
water tank hose and attach the hose leading to the drain 
in bilge area below the boiler. Then open the drain/fill 
valve and allow the water to gravity drain from the 
boiler. 


2. When the appropriate water level has been achieved, shut 
both the fill/drain valve and the vent valve on the 
auxiliary condenser. 


3. Energize the data acquisition system, computer, printer, 
and quartz thermometer power supply. Load the software 
program DRPML and check for proper operation. Then 
verify that the thermocouple outputs all correspond to 
the same ambient temperature. 
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Open the fill valve to the cooling water sump and adjust 
the flow rate to ensure the drain box does not overflow. 
The fill valve for the sump is located to the left of the 
boiler switchboard. 


Perform a leak test on the cooling water systems. 


a. Energize the cooling water supply pumps and adjust 
the rotameter flow rate from 20% to 80% while 
checking for leaks in the test section. Then secure 
the control valve and de-energize the cooling water 


supply pumps. 


b. Open the control valve for the auxillary condenser 
and adjust the flow rate to at least 30% and check 
for leaks in the system. After completing the leak 
check, reset the flow rate to the auxiliary 
condenser to the desired position, minimum of 10%. 


Place the system under a vacuum. To do this, shut the 
drain valve in the plexiglass cold trap and energize the 
vacuum pump. Check the pressure gage on the vacuum pump 
to ensure proper operation and open the suction valve 
located on the side of the auxiliary condenser. After 
the vacuum pump has reduced the system pressure to below 
3 psia, shut the suction valve and secure the vacuum 


pump. 


With the system under a vacuum, the heaters may be 
energized . Three switches must be placed in the "ON" 
position to energize the heaters. 


a. Switch #3 on panel p5, in the main hallway adjacent 
EOQ=soom H=k0 Ge 


b. The heater load bank circuit breaker on the left 
side of the boiler heater switchboard. 


c. The condensing rig boiler power switch, on the front 
of the boiler heater switchboard. 


When the heaters are energized, the power level should be 
adjusted to 50 Volts, using the silicon control rectifier 
on the boiler control switchboard. If the system is below 
2 pslia, the power level should initially set to 40 Volts 
to minimize vibrational shock to the system caused from 
vapor bubble formation and collapse. As the system warms 
up, the power can be increased, in increments of 10 
Volts, until the desired setting is reached. 
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After the system pressure rises above 4 psia, non- 
condensible gases need to be purged from the system. 
This can be accomplished by following the procedures 
outlined in step 6. Throughout the purging procedure, 
the horizontal tube coolant flow should remain secured, 
while the flow rate to the auxiliary condenser can be 
adjusted until all the gases have been evacuated. When 
the auxiliary condenser is warm to the touch everywhere, 
and condensate has formed in the plexiglass container, 
most of the non-condensible gases have been removed. The 
initial purging of the system may take between 15 and 30 
minutes. This process should be repeated every few hours 
during extended operation. 


To establish filmwise condensation on the tube being 
tested, the following steps should be performed: 


a. Allow the system vapor temperature to increase to 
=“90 C (at least 3800 microvolts on channel 40). 


b. Increase the flow rate in the auxiliary condenser 
to 506 or 60% and allow the vapor temperature to 
reduce to =75 C (approximately 3200 microvolts). 


c. Secure coolant flow to the auxiliary condenser and 
allow the vapor temperature to increase to about 85 
C (=3700 microvolts). This forms a steam blanket 
over the tube. 


dad. Initiate coolant flow of 80% in the single 
horizontal tube. 


e. Adjust the auxiliary condenser flow rate to maintain 
the desired temperature and pressure for data 
eg! 1 eerpion. 


If, after performing steps 9 a. through 9 e., some 
dropwise condensation persists, repeat step 9 again. If 
dropwise condensation continues, the tube should be 
removed, cleaned and retreated in accordance with the 
procedures listed in section IV. A. A cutout section in 
the insulation above the test section allows for a light 
to be shown on the upper surface of the tube in order to 
check for a continuous film on the tube surface. 


Run the software program DRPML by pressing the "RUN" 
key on the keyboard. The program will prompt you 
with questions for the necessary information 
required to run as listed by O’Keefe [Ref. 18] and 
in Appendix C. 


yy, 


Li 


Ze 


ales 


ae 


Ensure the system has been operating at steady-state 
conditions for at least 30 Minutes pricGmere 
continuing past the question "Enter the flowmeter 
reading". 


Vacuum runs are conducted at a heater setting of 90 


Volts and 1980 + 10 microvolts on channel 40. This 
corresponds’ to T,,, =48 €, and a Vapom velocTEy wea. 
Dem) Set 


The atmospheric runs are conducted at a heater 
setting of 175 Volts and 4280 + 10 microvolts on 
channel 40. This corresponds to Tco,. =100 Cyjauam 


vapor velocity of =1 m/s. Special care must be used 
when operating at atmospheric pressure to ensure 
Over-pressurization and rupture does not occur. 


When taking readings, double check the flowmeter 
reading prior to accepting any data point. If both 
vacuum and atmospheric pressure runs are to be 
conducted in the same day, the vacuum data should be 
collected first. This eliminates the long cool down 
time required after an atmospheric pressure run. 


After data has been collected, the system should be 


secured utilizing the following procedures: 


ie 


Secure power to the heating elements. Turn Ofbiieme 
Switches in the reverse order listed in step 7 of the 
Start-up procedures. 


Secure coolant flow to the auxiliary condenser. If the 
system is to remain at vacuum pressure until the next 
data run, then the auxiliary condenser can be used in 
assisting to cool the system down, provided the same tube 
will be used for the following test. 


Secure the coolant flow through the single horizontal 
tube by securing the coolant pumps. 


Secure the water flow to the coolant sump tank. 


To return the system to atmospheric pressure, slowly open 
the vent valve on the auxiliary condenser. To avoid 
contaminating the system, keep any foreign material away 
from the vicinity of the vent valve during this 
procedure. 


Ise an emergency should arise, such as an 
Over-pressurization or breakage, ensure the heater power 
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is secured first. Then, let the system cool down prior 
to checking for damage. 


C. DATA COLLECTION PROCEDURE 

After the tube has been installed and the system has been 
operating for a minimum of 30 minutes, data collection may 
begin. Each run consisted of data points taken at coolant 
flow rate settings from 80% to 20% and back up to 80%, in 10% 
m@erements. This allowed for each data point to be checked 
for repeatability during the run. Before a data point was 
accepted, sample readings were taken to ensure the coolant 
temperature difference, the overall heat transferred, and the 
Saturation temperature were stable. The temperature 
difference between the inlet and outlet coolant flow must be 
Semeilized (+0.01 C) prior to recording data. After setting 
the coolant flow rate, the aforementioned parameters required 


Berween 5 and 15 minutes to stabilize. 


D. DATA REDUCTION PROCEDURE 

For a Slice of the condenser tube, the overall thermal 
resistance from vapor to coolant can be represented by the sum 
g@mmemie COOlant side resistance (Ri), the fouling resistance 


Mrapeene Wall resistance (Rk), and the vapor side resistance 


ees shown in equation (4.1). 


Reocay ~ Kythet KR, +R, (4.1) 


a9 


Only clean tubes were tested so the fouling resistance was 


set to zero, therefore: 


R Serie nee (4.2) 


total 


Being convective in nature, the vapor and coolant side 


resistances can be expressed as follows: 


at 





Ro = ey (4.3) 
Be i ee 
R, — (4.4) 


The wall resistance term comes from the general radial 


conduction equation for a con@entrie cylinder 


_ 1nd,/D;) (aan 
3 Zines 


The total thermal resistance can be expressed in terms of the 


overall heat transfer coefficient (Ue ean the effective 
outside condensing area (A,), as follows: 
al 


Neosat ae UA, (4.6) 
Oo 


The effective outside condensing area of the tube can be 


calculated using the following equation: 


A, = ™D,L (4.7) 


Oo 
For the inside effective area, two different conditions 
are present. The area corresponding to the outside effective 


area, where heat transfer takes place radially, and the inlet 
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and outlet portions of the tube, where heat removal occurs 
meraliy. TO account for this “fin’ effect, the extended fin 
assumption and associated fin efficiencies were used in the 


following expression: 
Age = 1 Vie Bana bono (4.8) 


SUSSEITEUEIAG equattons (4.3), (44) weeatiGl “(46 )  -1nto 


equation (4.2) yields the following: 


E ae (4.9) 


Oye ys A; _ Hae 








From the data collected, the total heat transfer rate 
across the tube can be calculated directly from an energy 
balance using the coolant mass flow rate and the coolant 
temperature difference (from inlet to outlet) through the 


tube: 

Oehc ie Ty) (4.10) 
The total heat transfer rate can also be expressed in terms of 
the overall heat transfer coefficient, the effective outside 
condensing area, and the log mean temperature difference 
(LMTD) : 

QO = U,A,(LMTD) (4.11) 


where: 


(ao12) 
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The measurements for LMTD were obtained in the following 


Manner : 


® T, was measured from a quartz thermometer instaledea ea. 
inlet section of the coolant flow for the condenser tube. 


® T, was measured with a quartz thermometer I2msenee 
downstream of the mixing chamber at the outlet section of 


the coolant flow for the condenser tube. A correction 
factor was applied to T, in order to compensate. 
viscous heating of the coolant through the tube. The 


corrections, which depend on the inside geometry of the 
tube, the coolant flow velocity, and the type of insert, 
are listed in Appendix A. 


® T.., waS measured using the vapor thermocouple Siena 
40) 


Substituting equation (4.10) into equation (4.11) provides 
determination of the overall heat transfer coefficient 
directly from experimental measurements: 


yp a iteph 7) 


(4.13) 
oo aD) 


With A,, A;, R 


te and U, all known quantities, “Chemieues 


wi 
remaining unknowns in equation (4.9) are h, and h;, the 
outside and inside heat transfer coefficients. The modified 


Wilson plot technique was used to calculate these 


COect Eile Temuce 


E. MODIFIED WILSON PLOT TECHNIQUE 
In 1915, Wilson [Ref. 63], developed a method for 


indirectly determining the inside and outside thermal 
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resistance from an overall resistance. Wilson’s method 
requires the heat flux of the system to remain constant. 
Wome this application proves valid in theory, difficulties 
arise when trying to maintain a constant heat flux while 
varying the cooling water velocity, Marto [Ref. 64]. 

Khartabil et al. [Ref. 65] provide a comprehensive listing 
of modifications to the original technique developed by 
Wilson. For a detailed development of the modified Wilson 
plot technique, the reader is referred to Marto [Ref. 64] and 
Marto and Nunn [Ref. 66]. 

Using the modified Wilson plot method to determine the 
inside and outside heat transfer coefficients requires one to 
choose a general form for each correlation. The program DRPML 
enables the user to select from several options. The inside 
heat transfer coefficient correlation can be written in 


Simplified form as follows: 
jae GIN ®: (4.14) 


where 2 varies with the particular correlation used. Using 


the Sieder-Tate correlation [Ref. 49] (equation (2.19)): 


k 
OF = a Re* Preemie its oc (4.15) 
eb 
where the program allows the Reynolds number exponent (x) to 


be varied. Using the Sleicher-Rouse correlation [Ref. 50] 


fequation (2.21)): 
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K 
Q = sap + 0 OL Se rer ae) (4.16) 
where c and d are defined in equation (2.22) and (2.23) 


respectively. Using the Petukhov-Popov correlation [Ref. 51] 


(equation (2.24)): 


k 
© & Be) a Nee (4.17) 
D;| K, * K, (6/8) 17> (eee) 
where ¢€, K,, and K, are defined in equation (2°25) ee 
and (2.27) respectively. 
Seeley the outside heat transfer coefficient 
correlation may be written in simplified form: 
igen (A (4.18) 


o 


where Z can be selected from equation (2.1) (Nusselt’s theory 


[Ref. 23]): 


3 3 1/4 
nye 6) Oo ary (4.19) 
Cy a sige 18 





sat we) 


or from equation (2.62) (Fujii’s correlation (Ref? 25ia— 


z= st FS Keune (4.20) 


Oo 


F and Reo, are defined in equation (2.4) and (2955) 


respectively. 
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Substituting equation (4.14) and (4.18) 


(4.9) yields: 


1 - Fi Il 


= 7s 


R 
Oe a ce 











momen can be written as: 
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into equation 


(47215) 


(422) 


(4523) 


(4.24) 


(4.25) 


(4.26) 


Equation (4.22) can be written in simple linear form as: 


Y=mX + Db 


(4.27) 


Due to the temperature dependence of the parameters © and 


Z, the program utilizes an iterative procedure in the solution 


for each data point. A least-squares fit of equation (4.27) 


mieme@eceermines Cand a To calculate h,, the known value of 
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a@ can be used directly in equation (4.18); alternatively, with 
the known value of Ca) hh. canbe sea beuiac Goma sal equation 


ae 


(4.14), and applied to equation (4.9) in the following form: 


QD = 


.e) 








2 es \) (4.28) 
U, \ AA; a 
F. ENHANCEMENT RATIO 

The enhancement ratios in this thesis are used to express 
a comparison of enhanced tube results to smooth tube results, 
or to compare wire-wrapped corrugated tube results to plain 
corrugated tube results. 


From Nusselt’s theory, heat flux can be written as: 


ae (4.29) 


where: 








k a 2 h ie 
4. A ae ee asiog (4.30) 
H ¢ D, 
Heat flux can also be represented by: 
qe CN a Soe 


Equating equations (4.29) and (4.31) yields: 
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=e y “1/4 (4.32) 


sac WO 


Taking the ratio of the outside heat transfer coefficient 





Mmesmecanevenhnanced tube, (ho), ©O a smooth tube (h..), for the 
Swe i) 6 yrtelds: 
a a _ 
En = ee — aig = Pe (4.33) 
lee As a; 


Similarly, the enhancement ratio for a wire-wrapped corrugated 





Ble meOnstep ain acorrugated Cube ae the same (T.,--1,,./ Gan be 
represented by: 
A a a 
Eno = Ow? = WC = WC (4.34) 
lige Ao a. 


where the subscripts we and c refer to wire-wrapped corrugated 
and plain corrugated tubes respectively. 

For the purpose of this study all enhancement results for 
wire-wrapped corrugated tubes are based on €,, and wire- 


Deeeeea SMOOTH tube Enhancements are based on €,. 
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V. RESULTS AND DISCUSSION 


A. HEAT TRANSFER CORRELATIONS 

Lorenz et al. [Ref. 52] conducted a thorough investigation 
of several inside heat transfer correlations to determine 
which correlation most accurately predicted their data for 
turbulent flow of water in smooth tubes. Their findings 
showed that the Petukov-Popov and Sleicher-Rouse correlations 
(equations (2.24) and (2.21) respectively) were the most 
accurate in predicting the inside heat transfer coefficient 
over the range of Prandtl mumbers from 6.0 to 11.6. 
Consequently, O’Keefe [{Ref. 18] made a comparison of his 
results obtained using both correlations to process data, 
which was collected from the same experimental apparatus used 
in this study. His results showed good agreement between the 
outside coefficients obtained for smooth copper and smooth 
titanium tubes, when the two correlations were used with the 
modified Wilson plot data reduction technique. Based on his 
findings, O'Keefe [Ref. 18] processed most of his data using 


the Petukhov-Popov correlation. All data collected in this 


study were processed using the Petukov-Popov correlation so a 
direct comparison could be made with the results obtained by 


O'Keefe (Ref. 18 
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The data collection and reduction program, DRPML, listed 
in Appendix C, allows for the selection from two possible 
outside heat transfer correlations, the Nusselt theory or the 
Bigs Correlation (equations (2.1) and (2.6) respectively). 
Pore 5 Shows the reduced results for hj, obtained from using 
these two equations for a vacuum pressure run on a plain MHT 
Korodense tube. It can be seen that there is good agreement 
between the two equations. This indicates that the leading 
coefficient, a, obtained from the modified Wilson plot method 
using the Nusselt theory compensates for any vapor velocity 
effect inherent in the Fujii correlation. O’Keefe [Ref. 18] 
processed all of his data using the Nusselt correlation; so 
for comparison purposes, all further calculations of the 
outside heat transfer coefficient in the present work have 


also been made using the Nusselt theory. 


B. VERIFICATION AND REPEATABILITY OF RESULTS 

Several runs were made using the same tubes tested by 
previous researchers to verify system operation and data 
collection procedures. The five tubes that were tested 
included two smooth titanium tubes and two LPD Korodense 
titanium tubes tested by O’Keefe [Ref. 18]; of these, one of 
each was wrapped with a 0.5 mm diameter titanium wire. The 
fifth tube was a smooth copper tube wrapped with a 0.5 mm 
diameter titanium wire from Mitrou’s work [Ref. 14]. Figures 


6 and 7 are comparisons of the outside heat transfer 
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coefficient for a vacuum and an atmospheric pressure run 
respectively for data collected in this study and the study 
conducted by O'Keefe [Ref. 18] for the smooth titanium tube. 
The figures show good agreement for the two sets of data. The 
separation in the data sets seen in Figure 6 is due to a 
higher coolant water inlet temperature when O’Keefe conducted 
Wes runs. Tables II through VII contain a listing of all the 


data collected by O’Keefe, Mitrou and from the present study 


for comparison purposes in this’ study. The researcher 
initials are as follows: (0O)=0O'Keefe [Ref. 18], (L)=Long, and 
(M)=Mitrou [Ref. 14]. The average variation for a, the 


leading coefficient for the outside heat transfer equation, 
was less than 3.3%, with 90% of the results differing by less 
Ed@etie2>o. The C; values obtained in this study differed by as 
much as 8.5% compared with values obtained by previous 
researchers. This larger difference is possibly due to the 
uncertainty as to which HEATEX insert was used in O’Keefe’s 
experiments. Four different inserts were available and each 


had a slightly different length. 


C. ANALYSIS OF CORRUGATED TUBE RESULTS 

To establish a base line for the enhancement ratios of the 
wire-wrapped MHT Korodense tubes, an initial series of runs 
were conducted to determine the heat transfer characteristics 
of the plain MHT Korodense tube. The data obtained are 


presented in Table VIII. Figure 8 showS a near unity 


oye. 


Table II. SMOOTH TITANIUM TUBE AT ATMOSPHERIC PRESSURE 
















HEATEX Insert 


FONMAHT1 (0) 0.75 


FONMAHT2 (0) 










__FONMANT3 (0) 
FONMANT4 (0) | 













enhancement of the outside heat transfer coefficient for the 
MHT Korodense tube (having a groove pitch of 8 mm) over smooth 
titanium tube results. This is consistent with the findings 
of Rabas [Ref. 41] where, in a comparison of an extensive 
collection of data on spirally indented tubes he found: 

that there is no appreciable enhancement for 


condensation on roped tubes for groove pitches greater 
than 4 mm. 
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Table III. SMOOTH TITANIUM TUBE AT VACUUM PRESSURE 


HEATEX Insert 
Data Run/Researcher | ce | 


FONMVHT3 (O) 0.85 
O) 


FONMVHT4 ( 


0.79 
L) 


LSNTHV4 | 0.78 


FONMVNTS (0) 
| uswrva a) | 
| (L) 


LSNTNV4 





A comparison of the C; values for the no insert runs on 
the smooth tube (Tables II and III) and LPD tube (Table V) 
Igomicate increases in h, from 73 to 74% for vacuum and 
atmospheric pressure runs respectively. As discussed earlier, 
Wolverine [Ref. 2] reported a 25% increase in the inside heat 
transfer coefficient for the MHT Korodense tube when compared 


Paeie Pe Korodense tube. Companison of the C; values for the 


no insert runs on the MHT tube (Table VIII) and the LPD tube 
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Table IV. 0.5 mm WIRE-WRAPPED SMOOTH TITANIUM TUBE WITH 
HEATEX INSERT 


[Data Run/Researcher | CG; |e 
| rowmaner2 (0) | 2.39 | 99 
























(Table V) indicate increases in h; comparable to those 
reported by Wolverine [Ref. 2], 21% for the atmospheric 
pressure runs and 25% for the vacuum runs. This indicates 
that as the groove depth increases from the smooth tube (zero) 
to the LPD tube value (=0.33 mm) most of the increase in the 
inside heat transfer coefficient due to corrugation has 
already been realized. The additional groove depth of the MHT 
tube provides only a 21 to 25% further increase of the inside 
heat transfer above the LPD results (MHT tube groove depth 
=Q0.5mm), while the pressure drop doubles [Ref. 2]. 

Figures 9 and 10 show the effect of using the HEATEX 
insert on the overall heat transfer coefficient (U,)). The 
general trend for U, to increase with increasing coolant flow 


velocity is due to improved coolant mixing and a lower mean 
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Table V. PLAIN LPD KORODENSE TUBE 

























FONMANLT3 (O) 
LLNTNA2 (L) 


FONMVHLT1 (0) 
FONMVHLT2 (O) 







bulk fluid temperature. Increases range from 3% at a coolant 
flow rate of 3.83 m/s up to 9% at a flow rate of 1.02 m/s. 
This increase is attributed to the greater turbulence of the 
coolant flow due to the presence of the insert and can be seen 
mieeme higher C; values for the HEATEX runs for every tube 
tested. A comparison of O’Keefe’s data [Ref. 18] for smooth 


titanium and LPD Korodense titanium tubes, tested with and 


oy 


Table VI. 0O.5 mm WIRE-WRAPPED LPD KORODENSE TUBE WITH 
HEATEX INSERT 


Atmospheric Pressure 
Data Run/Researcher | CG | 


FONMAHL3T3 (0) 


(L) 
LLWOSTHAZ2 (L) 


| wonmvunst: (0) | 2.62 | 0.99 














LLWOSTHA1 











without a HEATEX insert indicates a 20% increase in U,) for the 
smooth tube and a 7.9% increase in U, for the LPD Cubemae 
tested at a coolant velocity of =2.5 m/s. For a coolant flow 
rate of =2.5 m/s the present study shows a 7.2% increase in Uy 
for the MHT tube when a HEATEX insert is used. This trend 
indicates that a more severe corrugation groove depth (zero 
for the smooth tube, 0.33 mm for the LPD tube and 0.5 mm for 
the MHT tube) results in a smaller increase in Up, with the 
addition of an insert. The grooves present in the MHT and LPD 
tubes cause greater turbulence in the coolant flow than the 
inside walls of the smooth tube. Therefore, the addition of 
an insert for a corrugated tube (which alxeady has turbulent 
flow) has a less pronounced effect on Uj) Compared §eomia- 


effect realized for a smooth tube. 


She 
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Table VII. 0.5 mm WIRE-WRAPPED SMOOTH COPPER TUBE WITH 
BeATEX INSERT 


pce mafeeeemee [qT 
, , 


FONMY H7ie2 
bow TCH va 
S71V296" 


“ Runs conducted with wire-wrapped insert 


























D. ANALYSIS OF WIRE-WRAPPED CORRUGATED TUBE RESULTS 

Seven different titanium wire diameters ranging from 0.1 
mm to 1.6 mm were used to wrap separate MHT Korodense tubes. 
Tests were conducted at atmospheric and vacuum conditions and 
the results were compared to the plain MHT Korodense tube data 
to determine the enhancement effects attributable to wire- 
wrapping. Wire-wrapped smooth tube data obtained by previous 
researchers were also used for comparison to establish if 
Similar trends exist. 

The enhancement due to wire-wrapping can be seen directly 
Miebigure 11 where the averaged U, (averaged for all data runs 


at each flow rate) are plotted for all the wire diameters. 
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Table VIII. PLAIN MHT KORODENSE. TUBE 










Atmospheric Pressure HEATEX Insert 
| patakun | cee eae 


LMNTNA7 


LNNTNAS 0. 80 


LMNTNAS 


Here the 0.4 mm wire-wrap appears to give the best heat 
transfer with an enhancement of about 3%. Similar trends for 
the outside heat transfer coefficient can also be seen in 
Figures 12 and 13. The plain MHT Korodense tube results are 
included to show the relative enhancements for each size of 
wire. Figure 12 shows an increasing trend with decrease in 


wire diameter starting from the 1.6 mm wire-wrap (P,,/D,=5) and 
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an enhancement ratio of 0.88 to the 0.4 mm wire-wrap 
(P,/D,=20) with an enhancement ratio of 1.08. Figure 13 
completes the data and shows the trend returning to the plain 
MHT tube results, Starting with the 0.4 mm wire-wrap and 
em@ang with the O0.lemm wire-wrap (P,,/D,,=80) . Two separate 
Figures were used to help eliminate the clutter from the 
overlapping data. feoke, PX seconeainomene saveraged, C wand a 
values for all the conditions tested (vacuum and atmospheric 
pressure, with and without the HEATEX insert) for the plain 
MHT tube and all the wire-wrapped MHT tubes. Ewos to -EOur 
tests were conducted at each condition to ensure repeatability 
of the data. 

The wire-wrapped MHT tube data were also checked to see if 
there was a relationship between the fractional coverage and 
the pitch to wire diameter ratio to the enhancement over the 
plain MHT tube (E,,.). Figure 14 shows the relationship 
between the enhancement ratio and the pitch to wire diameter 
ratio and Figure 15 shows enhancement ratio versus fractional 
coverage. Data obtained by O’Keefe [Ref. 18] for wire-wrapped 
LPD Korodense tubes are also included in the figures. 
However, the number of his tests are limited to only 3 wire 
diameters (0.5mm, 1.0mm and 1.6mm). The data clearly indicate 
that wrapping with a large diameter wire causes a degradation 
in performance. As the wire diameter is decreased, the 
enhancement passes through a maximum of about 8%, 


Perresponding €O a P./D,=20 and a fractional coverage of =5%. 
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Table IX. WIRE-WRAPPED AND PLAIN MHT KORODENSE TUBES 


Atmospheric Pressure 


— 


Plain 


3.27 


0.79 
| i.émm wire-wrap | 2.17 | 0.67 | 2.94 | 

| ) 

| 


0.3mm wire-wrap 






No Insert 










































O.2mm wire-wrap Zoe 


2.42 | 0.78 


As the wire diameter is further decreased, the enhancement 
ratio returns to unity, i.e. back to the plain tube 
performance. The results from O'Keefe [Ref. 18] for the LPD 
Korodense tube definitely show the same trends as the MHT 


tube, but the results are somewhat inconclusive due to lack of 
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data. His smallest wire diameter was only 0.5 mm. 

The enhancement results obtained for the wire-wrapped MHT 
Korodense tube were also compared to a wire-wrapped smooth 
tube. Marto et al. [Ref. 31] reported an optimum wire pitch 
to diameter ratio ranging from 5 to 7 for wire-wrapped copper 
tubes. Sethumadhavan and Rao [Ref. 34] realized a maximum 
enhancement at a fractional coverage of 21%, corresponding to 
a P,/D,=5. This value for P,/D, is much smaller than the 
optimum pitch to diameter ratio of 20 obtained in this study. 

To further understand this relationship, data were 
compiled for wire-wrapped smooth tubes studied at NPS and 
compared with the current data in Figures 16 and 17. Both 
sets of data were plotted with constant pitch curves for a 
better comparison with the constant pitch results of the 
corrugated tubes. It can be seen that the wire-wrapped smooth 
copper and wire-wrapped smooth titanium tubes with P,=4 mm 
have roughly the same slope. The same trend can be seen for 
the wire-wrapped smooth titanium and the wire-wrapped LPD and 
MHT tubes with a pitch of about 7 mm. Due to the limited data 
it 1s difficult to predict the optimal fractional coverage of 
the wire-wrapped smooth tubes for a constant wire pitch. 
However, it does appear to be less than the 21% obtained by 
Sethumadhavan and Rao [Ref. 34]. It can also be seen that the 
magnitude of the enhancements tend to vary with tube material 
and tube geometry. Mitrou’s data [Ref. 14] for wire-wrapped 


smooth copper tubes show results 45 to 55% higher than the 
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wire-wrapped smooth titanium tube results for an equal pitch 
(P,=4 mm) obtained by O’Keefe [Reft. 18) 2 This indicate saan. 
tubes made of materials with high thermal conductivities (such 
as copper) may provide additional enhancement characteristics. 
The effect of tube geometry can also be seen in the comparison 
of the wire-wrapped smooth titanium tube and the wire-wrapped 
LPD and MHT tubes (all with pitches =7 mm). The wire-wrapped 
smooth titanium tube achieved an enhancement =10% higher than 
the corrugated tubes of similar pitch. Figures 18 and 19 are 
identical to Figures 16 and 17 with the addition of the data 
obtained by Sethumadhavan and Rao [Ref. 34]. The tube 
material used in their tests is not known, but appears to have 
a thermal conductivity between copper and titanium (such as 
copper/nickel or aluminum). The most noticeable differences 
in their data and the data from research at NPS are the 
fractional coverage of the apparent maximum enhancement and 
the opposite slope of the constant pitch curves. This is most 
likely due the wire material used for wrapping their tubes. 
They used copper wire and all the NPS tubes were manufactured 
with titanium wire. The use of a highly conductive material 
may have produced more of a ‘fin’ effect with the wire and 
caused the peak enhancement to shift. Further study on the 
effect of the wire material is required to determine the 
effects. 

Possible effects on the inside heat transfer coefficient 


due to wire-wrapping both smooth and corrugated tubes were 
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Enhancement Ratio vs. 
Wire-Wrapped Korodense and Wire-Wrapped Smooth Tubes 


Figure 19. 


also investigated. Figure 20 shows a possible decreasing 
mee in Co with ancreasing Eractional coverage. However, the 
magnitude of the decrease is minimal and the overall variation 
fieeets withtm the uncertainty of the plain tube results. 
In summary, the enhancement of the outside heat transfer 
coefficient for a titanium MHT Korodense tube over a smooth 
titanium tube was close to unity. The primary enhancement in 
the overall heat transfer coefficient was found to be due to 
an increase in the inside heat transfer coefficient, 21 to 26% 
above the LPD Korodense tube coefficient and 110 to 117% above 
the smooth tube coefficient. The effect of wire-wrapping 
corrugated tubes provides minimal additional enhancement to 
the outside heat transfer coefficient, a maximum of =8% with 
a P,/D,=20 and a corresponding fractional coverage of 5%. 
H@wevyer, the addition of wire-wrap into the grooves of a 
corrugated tube may have greater benefit when used in a tube 


bundle arrangement due to improved condensate drainage. 
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ze) 
Fractional Coverage for Wire-Wrapped 


OF AE 
1 é 
Korodense and Wire-Wrapped Smooth Tubes at Atmospheric 


Figure 20. 
Pressure 


A. 


re 


VI. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


Enhancement of the outside heat transfer coefficient due 
to the corrugation of the MHT Korodense titanium tube is 
near unity when compared to smooth tube results. 


The corrugation of the MHT Korodense tube produces 
enhancements in the inside heat transfer coefficient 
From 110% to 117% above smooth tube results, and 21% to 
26% above LPD Korodense tube results. 


Wire-wrapping of MHT Korodense tubes yields minimal 
additional enhancement, =8%, when used in a single tube 
arrangement. 


The optimal fractional coverage for the wire-wrapping of 
MHT Korodense tubes is =5% (P,/D, =20) when used ina 
Single tube arrangement. 


The thermal conductivity of the tube material effects 
the enhancement characteristics of wire-wrapped tubes. 


Some relationship appears to exist between fractional 
coverage and the inside heat transfer coefficient. 


The modified Wilson plot data reduction technique 
compensates for vapor velocity effects when using the 
Nusselt theory in determining the outside heat transfer 
Soet ficient. 


RECOMMENDATIONS 


Manufacture and test additional wire-wrapped LPD 
Korodense tubes with smaller wire diameters ranging from 
O.1 mm to 0.4 mm to validate the trends found for the 
MHT Korodense tube. 


Manufacture and test additional constant pitch wire- 
wrapped copper smooth tubes with wire diameters greater 
than 1.6 mm and less than 0.5 mm and determine the 
optimal wire diameter for specific wire pitches. 


Test a set of wire-wrapped corrugated copper or 
copper/nickel tubes to determine the possible effect of 
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tube material on condensation, when compared to the 
titanium corrugated tubes. 


Develop a numerical model to approximate the outside 
heat transfer characteristics for corrugated and wire- 
wrapped corrugated tubes. 


Fabricate a set of corrugated tubes with outside 
electrical heater elements and conduct tests to validate 
the accuracy of the inside heat transfer correlations 
available for processing data on corrugated tubes with 
and without a HEATEX insert. 


Input the Wolverine inside heat transfer correlation to 
the existing program. Reprocess selected corrugation 
tube data and compare results with those obtained by the 
Petukhov- Popov correlation. 


Test wire-wrapped corrugated tubes ina tube bundle 
arrangement. Determine if additional overall bundle 
enhancement is possible with inundation effects 
included, and if the optimal fractional coverage shifts 
from that found for single tube tests. 


Fabricate and test a set of multistart corrugated tubes 
with and without wire-wrapping to determine if 
additional enhancements can be realized using the 
multistart corrugation geometry. 


Conduct additional tests on wire-wrapped smooth and 
corrugated tubes to determine if a relationship exists 
between fractional coverage and the inside heat transfer 
coefficient. 


Test the enhancement effects of various Wire Matenreee 
when used for wire-wrapping smooth and corrugated tubes. 


Investigate the possibility of bonding the wire-wrap to 


the tube surface and test the enhancement effects for 
smooth and corrugated tubes. 
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APPENDIX A. SYSTEM CORRECTIONS 


A. FRICTIONAL TEMPERATURE CORRECTIONS 

Coolant flowing through a tube experiences a temperature 
increase due to frictional heating. The magnitude of the 
temperature rise depends on the internal tube geometry and the 
velocity of the coolant. With the cooling water sump and the 
test section at ambient temperature, measurements were taken 
to determine the temperature increase due to frictional 
mmerecets at flow rates ranging form 10% to 100% on the 
rotameter. Two separate runs were conducted for the MHT 
Korodense tube, one with the HEATEX insert installed and the 
second with no insert. Figure A. 1 show the results for the 
two runs. The data was curve fitted to a third order 


polynomial as shown in Table A. 1. 
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Figure A.l Frictional Temperature Rise for the MHT 
Korodense Titanium Tube with and without a HEATEX Insert 
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Flow Velocity (m/s) 


Table A. 1 FRICTIONAL TEMPERATURE RISE EQUATIONS FOR MHT 


MORODENSE TUBE 
Inside Polynomial Equation 
Condition 
= 2.564x10°* + 6.263x1077V,, 


HEATEX 
Insert + 2.603x10°7(V ye + 7.830xX10 


W 
















No 
Insert 


6.) 790201 Cee wee 62051. 0 hrs 
PaCeA OSG) (Vv fo ee To tdi)” 


Ww 





= Temperature Rise (K) 


where: il 
V Fluid Velocity (m/s) 
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APPENDIX B. UNCERTAINTY ANALYSIS 


All derived quantities have an inherent uncertainty 
associated with them when their results are obtained from 
measured values. The uncertainty arises due to inaccuracies 
in the measurement of the physical quantities used in their 
calculations. The difference between the actual value and the 
measured value depends on the accuracy of the measuring 
device, the calibration of the device, and the experience of 
the operator. Eventhough the error associated with any single 
measurement may seem insignificant, it may have a profound 
effect on the accuracy of the derived quantity due to the 
calculations made. 

Kline and McClintock [Ref. 67] developed a method to 


estimate the uncertainty in the derived quantity, R, which is 











a function of several measuréd variables, x,, X>,))--- seen 
follows: 
W/Z 
OR pao aa dR, \’ (B.1) 
Ne = W. + | ——W. + W. ° 
(3E%) * [30%] cag 
Whee: 


Wp 1s the uncertainty of the desired dependent variable 


My, M5, +--+, KX, are ENG meastired sudependeneevantaem ss 
W,, Wo, ---, W, are the uncertainties in the measuhees 
variables 


84 


Georgiadis (Ref. 68] presents a complete discussion on the 
uncertainty analysis used for this experiment. Based on 
Georgiadis’ developments, Mitrou [Ref. 14] wrote a program for 
calculating the uncertainties in the mathematically derived 
Quantities obtained from data collected using the apparatus 
utilized in this research. Modifications were made to the 
program which allow for the data collected on MHT Korodense 
tubes to be processed. Sample uncertainty results for the 
following data collection conditions are presented in this 
appendix. 


i. Vacuum pressure with no insert and high flow rate (3.83 
m/s) . 


2. Vacuum pressure with no insert and low flow rate (1.03 
m/s) . 


3. Vacuum pressure with HEATEX insert and high flow rate 
(se e3-m/s) . 


4. Vacuum pressure with HEATEX insert and low flow rate 
mie 02 m/s) 


5. Atmospheric pressure with no insert and high flow rate 
i. 82 m/s 


6. Atmospheric pressure with no insert and low flow rate 
(6-02 m/s). 


7. Atmospheric pressure with HEATEX insert and high flow 
rate (3.82 m/s). 


8. Atmospheric pressure with HEATEX insert and low flow rate 
(02. m/s). 
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1. DATA FOR THE UNCER Tad tiie aie oes 


File Name: LMNTNU®8 
Pressure Condition: Vacuum 
Vapor Temperature = 48.574 Deguae 
Water Flow Pate (4%) = 68@.a@ 
Water Velocity = Bias m/s) 
Heat Flu« = | @6VE+O5 "W/m 
Tube-metal thermal conduc. = 21.9 (W/mK) 
Petukcv-Poopov constant = oT oe 
UNCERTAINTY ANALYSIS: 
VAR TABLE PERC ER Ieee Ra eeai lie 

Mass Flow Rate, Md 2.8 
Reynolds Number, Re he) 
Heat Flu«, q [eae 
Eso -Mean=len 02 pie Lae 
Wall Yesistance, Rw Oe sl 7 
Overaglrede | cae oS 248 
Wa teres teal leneenaaa Peau 
Vanor—-Sidegn. |. Cerio Gra 

2. DATA FOR THE UNCERTAINTY ANALYSIS: 
File Name: LMNTNVUS 
Pressure Condition: Vacuum 
Vapor Temperature = 48.691 (Deg se 
Water Flow Rate (%) = 29.00 
Water Velocity = pe 3 (m/s) 
Heat Flux = |. 14SE+OStaCw ma 
Tube-metal thermal conduc. = ra (ee, (W/o 
Pat ulov-P ogpwieeons tant = oe com 


UNCER GGIMIY “ANA yas. 
Ven PAB EE 


Mass Flow Rate, Md 
Reynolds Number, Re 
Heat Flux, q 
Log-Mean-Tem Diff, LMTOD 
Wall Resistance, Rw 
Overall. ].c 2s 
Water-Sideah. 20 se iat 
VaO0R = S102 SHi a) aCe 
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3, GATS FOR THE UNCEFTAINTY ANAL/SIS: 


File Name: LMNTHU8 
Fressure Comdi{$19n: Vecuunm 
Vapor Temperature = tees) 5 Deg | 
Water Flow Rate (4) = 320.09 
Wate ae Yt, = oes ied ee 
Heat Flu. 7: Pas25e705 Wm 2 
Tube-“etal tnermal conduc. = Ai Wim. kK? 
Peturer- onoyv cOnstamt = Sige Ud aa ee 
VHEEEIAINIY ANALYSI:: 

VARIABLE PERCENT SINGER TAINTY 
Wass = louo@=ate Md 0.80 
Feynolds Numnper, Fa 1.99 
edt. lis eG 1.33 
fog -leam—term Ort, LMTD ecte: 
Wall Pesistanca, Rw al sli 
Bera. Wane «0 Uo are 
Weer => eer.) eo OR Sis abe 
Vaeer swe mel. G. Ho carers 

4. DATA FOR THE UNCERTAINTY ANALYSIS: 

File Name: LMNTHVUS 
Pressure Condition: Vacuum 
Vapor Temperature = een Sys) (Deg C) 
Water Flow Rate (4%) = 22.22 
Water Velocity = Pave (m/s) 
Heat Flu- = eee OS aim 2) 
Tube-metal thermal conduc. = a ee (Wem.k ) 
Fethiowereney constant = Sel. 
PRESeTolNi yt earl slo. 

WieiniJtia)o' Gs PERCENT UNGER TAINTY 
Mass Flow rate, Md 3.0@ 
Reynolds Numoer, Re 5-203 
Heat Flu«, q 3.@8 


Beg -Mean-wen O1ff. EM1D ssi! 


Wall Resistance, Rw lst ei 
Overall ea. 1 26... Uo Soa ihe 
Wet Se or een) Stee. 1 Geo5 
Va or =o teen ie l.G «4 10 [ee a 


Sy 


5, OAlA For ThE CNet ANere Sp ele 


File Name: LMNTNAS 
PrassuresGenditian: Atmospneric 1@1 kPa? 
Vapor Temperature = Boreo od Ueoms 
Water Flow Rate (%) = 50.00 
Water el eiGet ty = ears @ 7 Sa 
Hea tee Ue = Jeet Oe Aan 
Tube-metal tnermal conduc. = es Wome 
Pet ov o00%— some lan = owes S3 
UNCERTAINTY ANALYSIS: 
alae) eles ble PERCEN PRUNGer lel: 
MasigacF beuiiea te. eric 0.80 
Reynolds Number, Re 1. td 
Heat hare ae 
Leqg-Mean= em Ut tama. sis 
Wall Fasistance, Fu eet 
Overall arenes... ave ees 
Wat 26 See eae cel «(ep mma Gar ais 
Vag ORs Sade ers lk ora aes Sid 
6. OATA FUR TRE UNCERTAINTY Sane sts: 
File Name: LMNTNA8 
Pressure Condition: Atmospheric (I@I kPa? 
Vapor Temperature = BFS, she] (Deg C) 
Water Flow Rate (%) = 22.00 
Water Veiocity = en (m/s? 
Heat Flus« = 5. IAS i 
Tube-metal thermal conduc. = Ra 9 (W/mK? 
PEevul sv Slo eee Olio amt = ..4653 
UNCERTAINT?’ ANALYSIS: 
VARIABLE PERCENT FUNGCERTALNT y 
Mass Flow Rate, Md 53.@2 
Reynolds Number, Re Sk) 
Heat Flux, q 5.04 
Log-Mean= fennel!) sei nell 
Wall Resistance, Rw oS / 
Overall eeueee eee S SWS 
Wate = Sige eee f Hi 8.04 
Vagcor-stge nel st = eo Pee): 
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7. DATA SOR TE LNCERTAINTY ANALYSIS: 


File Name: LMNTHAS 

Pressure Condition: Atmospheric ‘101 +Pa 

Vapor lemperature = adage + Deg C 
Water Flow Fate j %) = 88.00 

WSlern Siac tty - ao m/s) 
Heat Fiua = 4.Q@B8SEr0S WS te 
Tupe-metal thermal conduc. = ea 4) (W/mK | 
poe oy <r OD Opec Os Lamt = Ores 


DICE Tati y -ONAE YS 1&: 


VEBPIABLE PERCENT UNCERTAINTY 
Wass ir iew kate. Md 0.80 
Feynolds Number, Re Pere 
Hea tare oe. a 1.@9 
Lea=Mean-iem Oiff. LMTO Pod 
Wall Resistance, Rw ey 
eer la wt. 1 C+. 6 ec 
WAtCl =] shes. 1.6. Ht S205 
Vaeer-ovdewt 1 c. Ho 4.25 
8. DATA FOR THE UNCERTAINTY ANALYSIS: 
File Name: LMNTHAS8 
Pressure Condition: Atmospneric (1@0I1 kPa? 
Vapor Temperature = 99.9939 (Deg Cc? 
Water Flow Rate (2%) = 20.20 
Water Velocity = 1.02 (m/s) 
Heat Flu- = Seso(EtO0S ) CWem 2? 
Tube-metal thermal conduc. = ea ee (W/mK) 
Petulo7-Popey constant = Cee a Ae ee 
UNCERTAINTY ANALYSIS: 
VARIABLE PERCENT UNGER TANT y 
Mass Flow Rate, Md 5.@Q 
Reynolds Number, Re Cray 
Heat. fuse a 6.04 


Hegehean-rem Orff, LMI Pa ie 


Wall Resistance, Rw oe 
Overall H.T.C., Uo 5.05 
Watee eee.) Cet) ao le 
Wagors o.,ce co .1l.0 «~~ Ho Een gO 


ope) 


APPENDIX C. DRPML PROGRAM LISTING 


The data collection and reduction program, 


was used throughout this study follows: 


90 


DREMILF 


which 


hoo@e CS PML LONG, 

PoOey ss -EYISES FFOM ORPGr : rep ll Be Seer ee 

1Q0Q7 

Peer ) BE USED WITH MON-INSTRUMENTED TUBES ONL ¢ 

Byes) a8 —S GaAlTA Lh! THE FORMAT S& SWENSEN/O’KEEFE, LONG 
1010! CAN FEFFOCESS ONY NON-INSTRUMENTED CATA 

fOla! 

Borer is PROGRAM WAS USED fo COLLECT ALL THE NON- 

1014) INSTRUMENTED DATA TAKEN BY LONG (JAN-JUN 1993) FOR TITANIUM TUBES 
1015! 

rei7} MEANING GF ALL FLAGS IN PROGRAM 

1018! 

121d) ie ie Tee ares dels 

1Q20! (eye UF LO ee PR GSR AnM 

es |e! i le INPUT MGDE 

(ore Pee. | VALUE OF v1 USED 

oes! TFG: PONEO Oe SmogrTe 

Poe | INN: ise] ae 

3 | iWT: LOG NO. WITRIN PROGRAM 

1@26 | Beg) SeTESNATIVE CONDENSER TUBES 

12d! Ag: TSE. Oper raw 

1@23! ms: TUBE OLAMETER 

nee 3! Bee: PRESSURE CONDITION 

1030 i INF: OIMENSTONBESS Fille ReguirRed 

1031! IPE: PLOT F LRESREOULRED® 

POS. | IOV: OUTPUT REQUIRED 

ess | Cette PiSeer ate URRE LATION 

1d34! moe: CUTSTDE Hil. THEORY / CORRELATION 

PomoeeaM /Cc/ C(7) 

Boeoeeeonm 7CcSS/ T55(5) 

meee 7CcS5/ 15615) 

Bee econ yCcS// T5715? 

Pee OM *tcSe/ 158¢5) 

Peeve ooM “Fid/ [ft Istu 

meeeeolm Emf(c@) Twi) 

mweemeecoM ¢rr/ Yoal42),Timt42), Tfmr ,Ipe ,Qor 

Meee elie Wil/ Nrun,lim,luth,Imc ,I[fe,ljob ,Iwd,lfg,Ipco,lfte,Iwil ,Ihi ,loc,Inam,k 
cu,Re-p Fm Ax 

feel f/oeom/ OI Oe 01 ,Do,L,Li ,b2 

BeemeecOn ¢rric/ I[stuc,Inn,Iwtho ,Vw ,Vws ,I namo 

meee orn ©.1008609) 25/27 .94469 ,-767545.8295 , 75025595. 81 
Reem —-3247486569 €.97666F+!1 ,-~.6519ZE+13 ,3.94078E+14 
fees EAD Cts ) 

Been 2fo.15,2.5943e-. =F 26 71E-7 ,3.c941E-11 ,-9.7719E-16,9.7121E-20 
base READ T5S5( +) 

Pemeebete o1o.15,..5878E-2 -S.9853E-7 -3.1242E-11 .1.327SE-14 ,-1.0188E-18 


Sal 


PEs e See 


DATA 275.15 ,2.S923E-2 ,-7.39353E-7,2. 86255-1012 S7 7 ep re 


PEeR aise ase 


DATACL 13. 1S) ee Se eo 


READ TS8(*) 
or]. vlo5 to | 


Ossp=.'524 | 


be =. 054325 | 
mewUniehctec Is o3) 
Bio 


PRINT USING “4x, 
Pe (US [NG ea ae 
PRIM) We NG ae 


INPUT Iss 


rise) 


THEN 3094 


Inside 


? 


-o7-E=/ ,4,.Q657E-11 =1.2791E=-15 6.400 eee 


Outside drameter sf she cutter ean 


Giameter of stainless steel test section 


en 


d 


ih) 


O 


it 


7 
Ad 


Condensing length 
Inlet 
Outlet end 


Fin length" 
“fin length" 


Gt cn os 6 48 
Cata or re-process previous data’ ” 
N Analysys' “" 


TIMPUT “SELECT FLUID <2=WATER | =h=113. 5 —=6 yea 


Nie ire dee 
Slelale 
Tjob=@ 


INPUT “ENTER INPUT MODE “@=S0540 > =F Tees oar. 


Im=ime | 
S/S|e/< 
IF Im=t THEN 


INPUT “ENTER MONTH, 
QUTFUT 729; "10 =daces 


BU RReDe SeniSE aa 
ENTER 709;Dated 


ENO 
LP Ul jcb= (ane 
BEET 


DATE AND TIME «MM: O02 Hee Wii Sse ee oer 


INPUT “Sf IP PAGE “ANS Hl 7 Bilbo raael: 


ENO IF 
PRINTERS DS ul 
TF Im=t THEN 


ENTER 709;DateS 


eAnUyN ap 
= N= 
PRINT 


PrN US tNG eae 
[hos 


IF [jyob=t THEN 


"NOTE: 


Month, 


P 


date and time :";DateS 


rogram name OR aes 


oa 


O35 Cm ~} 


te 


GS - 


folie ¢ 3 6) fort sp 
SL a eS 
S10 (sor rs — 


f 
L 


Beer 


ENPUT 


hie 


[m=] 
EEEP 
LNPUT 


“GIVE A NAME FOS THE RAW DATA FILE” 0 files 


SELECT 4 
THEN 


PRINT USING 
PeeateE GOAT O files 52 


feetoN SFite 16 O-filss 


SEER 
INPUT 
Tnn=0 


TENIEP GEOMETRY. COBE (r=FINNEO G=FEAIN® {lia 


Pebiter 15-1 


Beer 

PRINT 
PRINT 
BR INT 
ee INT 
eT 
INPUT 


Inn 


ENE eee m Wee” 
BOEr UE is 


CieaO=r IND }=S1ORED Mee ea 


“16% 


@=NONE 
l=TWISTED 
C=HWIRE 


ZsHEATEX" 


Wee oT CGFi let 1 te..inn 
FOUR UNINSTRUMENTED TUBE 


Iwt=0 
Fh=0 
Fpo=0 
Fw=0 


Istu=0 


Istuo=[stu 


IF Ifg=0 THEN 
"FIN PITCH, HEIGHT ANO WIOTH IN MM, Fo Fh.Fw" Fp Fh Fu 


INPUT 


pe | 


"File name 


al gl am 


WRAP “ 


ee uy SFrle:ilwt Fo ,Fw,Fh 


Bok 


BEEP 


peo NIER > | 


PPINT 
EN | 
PRONT 


SRUGENT & IDATA TOsbER Ree e0CEs se = 
V-=SJENsEN7 0 WEEEE, LONG <DERPAVET 
PVN Pet beN7 Mit ROU /COUMES7GUTTENDORE © 


INPUT Istu 
Istuo=Istu 


Beer 
BRN 


IF Istu=1 THEN 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


Boe 


S LUDER NAME: ~ 
Q=VAN PETTEN” 


~ 
a 


=MITROU" 
01.0181 | =e 


3=GUTTENDORF " 


35 


oleae Led 


PRINT USING “!5X% °°" This emnalycsis Gomme oneni1 2 


OF THE 


USING 


USING 


USL NG: 


BS TING 


USING 


USING 


PRINT ~=SWENSEW 
PRINT S=O'KEEFE" 
PRINT 6 =LONG" 
ENDO IF 
INFUT inam 
Inamo=inam 
BEEP 
INEUT GIVE TRE va e 
PRINTER [S 791 
IF Inam=@ THEN FRINT 
IF Inam=!t THEN PRINT 
IF Inam=2 THEN PRINT 
IF Inam=3 THEN PRINT 
IF Inam=4 THEN £&-iNT 
IF Inam=S ‘HEN PRINT 
IF Inam=6 THEN PRINT 
PRINTER IS 1 
Beer 
INPUT 


ASSIGN @File TO D_ files 
ENTER: @Pilect foe ian 
IF Istu=@ THEN 


ENTER ACP ike lit ee eureeh 


Giese 
IP lfe=0- THEN ENTER SCR ulema 


lie 


Bibel F 


ENGIIE 
IF Ijob=1 
IF Ift.@ THEN 


GEEE 


PR UNRGR iS 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


lwtho= 


USING 
WSING 
USING 
USING 
USING 
[wth 

[wth 


lfg=l 


THEN 


foo 
1345 


EXISTING SATA FICE Otis 


Bees 


NS 


a Vio), eran 


Maan 


CPS) sae 


"Data 


* “Dia ta 


“Data 


“Data 


Data 


Data 


""Data 


THEN ENTER @File;Fp,Fw,Fh 


"ak, “Select tube type: = 
"6X ,""@ Copper Tube""" 


Bahra a 


"6% = SMOOth iitanium. | ube 


"6X os 


94 


taken 


taken 


taken 


taren 


taken 


taken 


taken 


"ENTER NUMBER OF DATA SETS ss (OR Rare rien 


Dy 
Dy 
by 
Dy 
by 
Dy 


by 


Wolverine Korodense LPO Titanium Tube" 


it 


"" 100A “;,.005% es 


"3 Wolverine Korodense MHT Titanium Tube""" 


VAN PET ig 
MI TROO 


COUMES” @ 
GUTTENDO| 
SWENSEN* 
O° KEG ries 


LUNG” “@ 


15.26 
Sa re 


oe 
ens! 
ie 
13350 
Pos 
ee 
aos 
1334 
Lo5c 
S65 
hogy 
oc 
Koei 
1348 
1341 
[2 .. 
1343 
1345 
1348 
1349 
eS 
iss. 
(3a) 
1360 
463 
Stags 
to67 
1363 
Lage 
ete 

Load 
lait 
1384 
les9 
14Qc 
1405 
1457 
1458 
1459 
1468 
1461 
1462 
1463 
1464 
1466 


BEEP 

PRINT USING ‘4% ''Serect tube Enhancement used: """ 
pein USING 5. ) -MOOTH Tuce: 

eee eto aCe es FINED [WBE ” 

PRINT USING “64° '2 wIRE-wFAPPED SMOOTH TUBE"”” 

PRINT USING “64, °"2 LPD FORODENSE TUBE’ ’” 

peti estING “6% “4 WIFE-WRAPPED LPD KORODENSE TUBE*”” 
Peegva west NG "6x ~°S MAT rOPGDENSE TUBE ©” 

PRINT USING "6X ,""§ WIRE-WRAPPED MHT KORODENSE TUBE"”" 
INPUT Ityp 

PRINTER [5 721 

BEEP 


PRINTER [IS | 
PRINT USING “4x 
ee int USING “bd 
PRINT 
PRINT USING 
INPUT Ime 
PRINTER IS 1 
BEEP 

Itds=| 

IF Iwth=@ THEN 


USING "6x," 
yy area caine 


PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


USING 
USING 
US LNG 
USING 
Itds 


eno iF 

BeeNicR [5 791 

IF Iwth=@ THEN 
Or. 0107 
Do=.019@5 

END IF 

IF Iwth=1 THEN 
Di=.01347 
Do=.01585 

ENG IF 

IF Iwth=2 THEN 
Do=.@1585 
Di=.@1386 

Eno oir 

IF Iwth=3 THEN 
Do=.016@7 
Di=.@1355 

ENios ie 

01=.01905 


"Select Material Code: ” 


"@ Copper 1 Stainless steel"*” 
ey ome Aum OF sole. Cu NE 
Titanium 


"4X ,""SELECT TUBE DIA TYPE:""" 


ea. 2 SmMAb 
foe  ) “MERU CER AUET I~ 
Cae VEG 

YO OF MEDIC ata ek (1 UBES 


POOP Ore MEO TUM “Upc 


25 


1467 
1468 
(eusys, 
val 
1474 
L477 
1478 
1479 
1484 
1485 
1487 
1489 
eon 
1495 
nel 
1498 
1501 
15@4 
fed 7 
'SO8 
foi 
S57 
1538 
15453 
1S46 
iS49 
[SS2 
| S56 
(Sar? 
ifsc 
(253 
1S60 
S60 
Heileye 
iSoc 
1564 
eos 
1567 
f5o6 
ieSo.g 
1S7@ 
eal 
Se 
owe 
1574 
[S¢S 


Do=,.01557 
IF [tds=0 THEN 
o= soe? 
Ji =s003Se5 
ape. le 
if btds=c Ve Go) ges 
IF twth=1 GR iWtn=5 THEN oS se 
IF Twth=| OR ITwth=3 PRHEN G2-—.21ee5 
iF [wthe=2 THEN G61=.015S37 
IF Iwth=2 THEN O2=.01567 
IF Imc=Q@ THEN Fou=385 
IF Imc=l THEN Fcu=16 
IF Imc=2 THEN *cu=167 
IF Imc=3 THEN Feu=45 
IF Imce=4 THEN kecu=l!.2 
Rm=Do*EUGt Do Ur C2 <wer Wall resistance based on outside area 
Seer 
INPUT "ENTER GES SsUee COND oh Yay =a Sloe 
laeo= lac 
Inft= 
BEER 
[fe=l 
PRINTER® iS 7C! 
PRINT USING “l6xX,“"This analys1s inetmces end— tame ce meee 
PRINT USING “16% ,"“"“Thermal conductivity = "" 30.0,"" (W/m lk eee 
PRINT USING “l6K ,""Imside diameter, D1 = "" 00.00,"" (mm )” “Sgn 
PRINT USING “16x ,""“Qutside diameter, Do = "" 900.0D,"" (mm)> sg 
Thi=0 
epee sisi 9 is 
PRINT SELEG? INSIDE Cee elay ironies 
FRINT J=S1ESER-—lA ee Obes eee 
PPINT 1 =SLEICHER-ROUSE- 
ee Tha S=FE IU RGvV-- Orel 
Cie ie act 
IF [hi=@ THEN 
Beer 
INPUT “ SELECT REYNOEGS, EXPONENT Speco 
END IF 
Toc=0 
Gee P 
PRINT 
PRINT SELECT OUTSIDE THEORY/CORRELATION FOR WILSON ANALYSIS: °” 
PRINT G=NUSSEET TREORY= (Geet aie 
rile \=FUUT? (Sage CORREA nol 
INPUT loc 
BeEe 
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IS’ = [tm=l 

foe FRINT 

Mego) FPINT Sewee Se eClaNnl TEnrPerAlene Se MeesurEMENT >” 

Poem [situ=0 THEN PRINT G=S0 (GV repeat 1. C" 

1SeQ PRINT V=(UAr Te VHEeMOMETER, “DERavET 

Pssl PINT 2=10- (UNC LION THERMOP IEE: 

| aysem NPUT Itm 

Peoom FRINTER IS 7O1 

eee Ltm=9 THEN PRINT USING “J5X,"“"This analysis uses the SINGLE TEFLON T'C 
Peacings’ " “ 

Peele ttm=! THEN PRINT YSING "16X,.""This analysis uses the QUARTZ THERMOMETER 
readings" '" 

Reece itm=2 THEN FPARINT USING "I164%,""This analysis uses the !@-JUNCTION THERMO 
eieeereadings’”” 

Peo 7eeiic=! |! FOR MODIFIED WILSON 

Bogemeie 1hi=8 THEN Ci=.007 

Meee [hi=!l THEN Cie=l. 

eee TP hi=2 THEN Ciel. 

Moogee.r [wil=1 THEN 

1601 Seer 

1602 PinoieeeNTER @: IF OLFFERENT FROM STGRED VALUE” ,Ci 

hegee END IF 

1604 PRINTER IS 791 

Poe 1[hi-O THEN PRINT USING “16X,""“Modified Sieder-Tate coefficient saree 
me ul 

1606 IF Ihi=@ THEN PPINT USING “16% ,""Chosen Reynolds No. exponent = a) 
DD" ;Reap 

15@7 IF Ihi=1 THEN FRINT USING “16X ,“"Modified Sleicher-Rouse coefficient = 
ee 301 

Pegeeeerr ihi=2 [THEN FRINT USING “16X,""Modified Petukhov-Popov coefficient = 
ree oo! 3 C1 

1609 fF Inmn=@ THEN PRINT USING "16X%,""Using no insert inside tube""” 

poet emm=2 THEN FRINT USING “16X,""Using wire wrap insert inside tube"’™ 
eemeetee ian—> THEN PRINT USING “16X ,""Using HEATEX insert inside tube™"” 
Peybeeeir [stu=2 THEN 

reas Patna= WHEN PRINT USING “16X.°°Using twisted tape insert inside tube™ 
eo.) ELSE 

els eens) ToeEN PRINT USING “~16X °°Usinmg wire wrap insert inside tube” 
feos) END IF 

1619 IF Iic=@ AND [fe=1 THEN Ac=26.4 

PaeveeelF fic=!) THEN Ac=Q. 

ices BEEP 

eeoe lF Lycb=| THEN 1648 

feeon PRINTER IS 1 

eee INPUT “NAME FOR A FEMPORARY PLOT FILE (TO BE PURGED?)” ,P_files 


Sy 


[625 Pltilet= Curry 

[B25= "SEER 

1634 CPREATeG BOAT Perera so 

Vot2 S515 SFriee eo ee wees 

bede OP Vien =i" Tae 

Les lo.= 

1654 G0TO wi 6s 7 

[654-2 ENE ae 

i6ee Sees 

1661 INPUT *~SEREST SGU PUTS 2-=SHORT. 1-20 G. aeeee. 

(Eba. Vbov=tev-] 

1667 PRINTER [5 701 
672 IF It,p=0 THEN PRINT USING “l6A4 .° tube@eanancenent SMOOTH 1Ue= 

(675 IF Ityo=!| TREN PRINT USING “I6A%.° Tune Epnaneenent FINNED ?06e 

[o7A [Et ye.) THENG EE ES UNG Se. Tube Enhancement WIRE-WRAPPED sMOOT re 
TUBE" | 
l6°S IF [tyop=s5 THEN FRINI USING “164, Tube Enhancemen: LPO +ORODENSE Jeo 
16°58 LF Ityo=4 THEN PRINT USING “I6X,""“Tube Enhancement WIRE-WRAPFED LFO KO 

RUCENSEIUEE 

1679 TF Ityo=5S "REN PRINT USING lex “Tuge eanamecenen MHT KORODENSE SS Veae 
166@ IF Ityp=6 THEN PRINT USING "I6X,""“Tube Enhancement WIRE-WRAPPED MAT K@ 

EOCENE sucess 

(S62. 95e Er 

1633 IF Imc=@ THEN PRINT USING “16X."*Tube material COPPER 

1664 IF Imc=| THENSERINT Wists ie ee lube tmaen sat STAINLES S= oS ianee 
1l6G5 iF ime=2 THEN PRINT USING “1l6A> (time wmat er wa ALUMINUM as 

1666 ir [Imc=5 THEN PRINT USING “16X," "Tube material 90/10 CU/ hi 

16? IF Imc=4 THEN?PRING USING Sex, Tube mater ral TITANTU 

i688 IF I[pe=0 THEN. PRINT SUSING | 16k Presstemescnce 2m VACULT aa 

1669 IF Ipc=l THEN FRINT USING “16X,"““Pressure condition AIMOSPRER la 

159Q0' FRINT USING “16% ,""Fin pitch, width, and height imm): °” DD@SD 32 eee 
OR Sle) Jan ala 

1691 IF “\lwil=0 UR [Waske2) AND Im=2 TREN 

632 Tjob=! 

bas 2 Iwd= 

1696 Cab, Wal someG 

1633 ENO IF 

170206 S =O 

vc Slee = |S ten 

eee PRINT 

a ae US Waite Miele 

1724 PRINT USING 10x" Davaue. Uo Ho Qp et 

Ts Rexp""" 
lkZeS PRINT USING “1Gx °° teCm7 5 Cie ee 7 ei (W/m ae ae 
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ir, hit = 
= CO ime) vr 


ye OSL 400) Teo) To" 


| toieear 


nN 


(=) 5 [i (Wy Wis TR RON fo OS) & a fs) fe (it 


ma il oO 


ELSE 
PRENIWERiG See '"Cata was Us 
cians 
St as [ys OX ¢ m/ 5 Nam W/ 


m=! THEN 
Beer 
PePUP EIKE TO CHECt NG GONCENTRATION (1=Y .@=N)?7" Ng 
IF J=t THEN 
MAVPUT 709; AR AF4AG AL41 URS" 
CUrree! "9S; (AS 5A” 
END IF 
Beer 
INPUT “ENTER FLOWMETER READING” Fm 
OUTPUT 7093; "AR AFGO ALE2 VRS" 
Meet 709s As SA" 
ENTER 709;Etp 
meeeuyT 789: "AS SA” 
Beoe 
Pew Ph  EONNECT VOLTAGE LINE” ,Ob 
Bek 7J9:6vol 
BEEP 
PePuy ~OPSCONNECT VOLTAGE LINE” .Ok 
IF Bvol<.1 THEN 
BEEP 
BEEP 
iv eV Tee VOETAGE, TRY AGAIN!” Ok 
Erne, asia 
eet F 
Pe PUl Oss AS. SA 
ENTER 709;Bamp 
EHO. ie Eto 


oo 


DUTPUT FOS eee hes nes ae vine 


1S65 

1374 Nn=? 

395 “OR $=9 TO ian 

1g73 QUTFUT T@39:"SS SA 

165 a= 

1583 FOR t=! TO 10 

1@9l ENTER 7@9;E 

1394 Se=setE 

ens ie ies has 

1980 Emi ( | =Abs( Se71Ge 

1916 Emf( I )=Emf ct )*1.€+6 

Jashies: NEXT I 

ba COMET eas seamen 

ged OUTPUT 7152) ieee 

aes He | 

Va20 ENTER. fist 

mes QUT al oa wee ee 

| Sisis Wel we 

es) Shs ENTER 713372 

base CUP oils Thee 

{94S Napili: 

1948 ENGR eaeeniee 

Peas TAS ite 

Sis! OUT 715 ibon oe. 

196@ BEEP 

| see INPUT “ENTER PRESSURE S@GnGe Reo O DIG SU rganeeees 
1971 Pvapl=Pa@a*be34- 72) Psi seera 

boa OUTPUT 709;"AR AFE4 ALG4 YRS” 

sre OUTRO a7 OS Seon | PRESSURES TRANS@UCER 
ba74 S$s=2 

1975 FOR K=1 TO 290 

[ses ENTER OS. ee eeon 

S77 Ss=SstEtran 

sifaei NEAT K 

aa Ptran=ABS(Ss/2@) 

1980 Beer 

sisi PRESSURE IN Pa FROM TRANSOUCER 

[shes Pyap. = (-f2.55b04*e tran 14.757 Soe 
| Setse Slee le 

1986 IF Istu=@ THEN 

spans ENTER @File;6vol ,bamp Eto -m il 2  -vaol Pyveane sere ae 
1990 Eis e 

poae ENTER @File;Bvol ,Bamp ,Vtran,Etp ,Emf(@) ,Emf (i) Emft(2) ,Emf(3) emf (4) ,F 
mi | 1 Te ePhoeeiater 

sisi! END IF 

1356 IF J=l OR J=2Q0 OR J=Nrun THEN 


5810/8) 


Posy Ng=! 


bags BESE 
| Sisk Ng=0 
~800 S|) Olas 


meee END IF 
Beeee IF Istu=0 THKEN 


2.008 Tsteami=FNT.svS7iEmri@)) 
~@@9 Psteaml=Isteam!-2735.15 
2012 Tsteami=FNTvsvS6(Emril>) 
ea isteame=Isteamei-2,3.15 
EOE Tsteam=lsteam! 

eos Troom=FNIvsvS8¢ Emr (2?) 

a ae feaom=!rocm— 7 ..."S 

2238 Meonm=FNivesevos(Emft7 i) 
<@239 Heom='con=.:4.'5 

moe.) ELSE 

E04 5 Peteam=PNivsvl Emit 07: 
-O44 Pecem=eNivsv  Emte 2 )) 
Bvas Teem-PNivsyvtEnri 42) 


moss END IF 

2048 Fsat=FNPvst(Tsteam) 

meer ong—| 5529-122*( Troom-26.85 )/50 
2053 Rowater=FNRhow( Troom) 

-@63 IF Istu=0 THEN 


208 | Ptestl=Pvapl 

BOB. precti=rVapc 

moeo ELSE 

2O84 Ftestc=(Phg*Rohg-Pwater+Rowater )#*9.81/1000 


eos END IF 

2087 Phs=Psatel.E-3 
mOeewerho=-rtestcel .E-3 

Bvege fl t=Pkis 

“O92 STsat=FNTvsp( Psat) 

2098 VYVst=FNUVsttTsteam) 
Beoeeeerpno=(Ptestc-Psat )/Ptest2 

mee Past =!-Ppng 

eee Mwv=!18.016 

eleS IF Ift=1 THEN Mwv=137 ir BE "CORRECTED 
Oreo Lft=2 THEN Mwv=62 

bese Uing=(Ptest2-Psat )/Ptest2 

Seer Tenor! /( 14 1/UFng-| d*Mwyvs28 2.97) 
clo? Mfng=MFfng*1@2 

2-128 Vfng=VFng+12 


Meee) BEEP 
mimeete IF lov=c THEN 
Enon PRINT 


Jeo 


zee 8 RECORG TIME OF Tara sGeoate 

es IF Im=1 THEN 

2142 CuTPPS 20S. ae 

Lidl ENTER SOs leds 

fe ENE 

144 PRINT USING “iOxk ,°"“Data set number = “" 00D ,4%. 1443] ee 
BUS OUTPUT 309%) AR Aes ee URS. 

Seo OUTED Sse on 


Ses) Np) ss 
2152 IF lov=2 AND Ng=1 THEN 





gels PRINT USING. ove msi Ptran Tmeas heat NG %""" 

Ren Pals: PRINT USING “10x ,"" (kPa? (KPa) (C) Ge Molalla 

eo PRINT USING “114,1¢20.0D,047,1(20.00,5%) ,2( 30.00, 2X ) 2X , | ( MSD Se 
Pkp ,Tsteam Isat ,Mfng 

2164 Peulue 


2167 END IF 
217Q IF Mfng .S THEN 


Lar ae Beer 

tgs [F Im=! AND Ng=l THEN 

pea PS Beis 

o18e PRINT 

@ els PRINT USING “10x ,""“Energize the vacuum system “"" 
Blog BEER 

ee INPUT “OR PO ACCEPT THIS RUN SO1=¥ Dh) ee 

cule [IF Gk=0 THEN 

eee, Beer 

. 280 BISP “NOTE sels DATA SSET Wile Ge) Seana 
2cOS Ween 5 

UA S GOTO 1782 

megs eNOS 

re ENG TF 


~-1S5 END IF 
~2lt& IF Im=! THEN 


ao IF Fm.t10 GR Fm 100 THEN 
a ee: Ifm=Q 

ui hele es 

a [INPUT “INCORRECT Bie leO0 Ger ls Qatar e = erm 
CAG 5 IF I fm=@ THEN 1804 
2256 ENG te 

2299 SEND ie 

2242! ANALYSIS BEGINS 

~2-43 IF Istu=@ THEN 

false TiLl=FNTIvsvS8(Emf(3)) 

EN ae Tirc=FNTIvsvS5(Emf(S) ) 

Fao Tol=FNTvsvS8(Emf (4) ) 
ease Tol=FNIvsv55(Emf(6 )) 


LOZ 


25S al ipa diene eer. 155 

20. Ga, eg ea Ri 

oye iol = eames oS. 45 

ie. ie =) Oo ee, SS 

NS Bae = fo] oti 

ea C eo 852 lie 

SoS Tdelt=ic-T] 

R255 Etpl=emr(3)+Etp/2Q. 

bos4 ee oe 1b -2 |  SUdbIE-S*Etel tt cl FO IE=-10*Etp1 2-5. 164E-15*Etpol 2+23. 
Poe | OF ip| 4 

foo watoace Cete.o/\0. 

eo a3 Poo= it i Ple 1s 

feo fei av=2) 1HEN 

2 5: Pe ies NG. ex OTN TOUT] TING Tees Cee BELTS 
Elle 

eet Boles NG it - JFEFLON ) CHOAR Ie 
eo64 Beet NO OX eo UD 24) 21 rl tol lle stdeli ldel3 tris 
Loos =o es 

Soo 7 Er i=ABS(T11-TI ) 

easy e Beocroott1.-1 1) 

A aieks, pranterm 15 | 

2570) BEEP 

e575 Bee =nbsit P2-)1)=( Tris) )/( 12-11) 

Oe eee OS eaNee im= | HEN 

Eegec! BEEP 

ae eee) . QC) ene TeP ite CIFPRee BY TORE THAN 52° 

~ 380 Oke=| 

ao | Heol e=0 One Ee. .25 AND m=) THEN 173@ 

Boo c END IF 

2565 PRINTER [5 791 

eae ELSE 

Tess S peter NT ysevCEmt (2 3 

B55 Seac-rNGrad’ «1 il+Tlc7*.5) 

Boo / To=TitABS(Etp)/(10*Grad)*!1.E+6 

2339 Prv=Ty 

boo S To3=To 


Peg eNU IF 

e239 6CdTF: GUL stu=@ AND Itm=0 THEN 
2O99 Tlizftil 

mood Teo=Tol 

Sege. END IF 

ooo iF I[tm=! THEN 

Bead elete= ia 

2540 oe 12 

ese END IF 

2200 IF [tm=2 THEN 


OS 


——4 


EN Ore 
lav Q=0 Tile eee 
[ft=d 
Cow=FNCpwit Tavg? 
Rhow=FNRKow( Tavg? 
IF Istu=@0 THEN 
Md= (S27 40Se moe. UO. 
Mda=Md*(1.0365-1. S6644E-35<] 115-2526 -6 ee eee 
spe |= 
Md= 1) 204 60SE—276 2609525 em 
Ma=Md*C IW565-1. sbe44b—-Se 55252656) eee eee 
Spares 
Mf=Md/Rhow 
QO ita ee ee 
Jws=Vwe(Oi/i|. 276-2?) 2 
VeGimee i i ete cee! 
1.252 eer 
N=Ma*eCows | eco- tli) 
Qp=Q/(PI*DoFtL) 
If t=Q 
KW=FNKW( Tavg ) 
Muw=FNMuw( Tavg) 
Rei=RhowtVwtDi/Muw | ASSUMED SAME FOR KORODENSE 
Prw=FNPrw( Tavg) 
Fel=Q. 
Fel=Q. 
Cf=1. 
Preuf=Prw 
Fer:f=Re} 
iF Ini=@ THEN 
Ome=Fe1 Pe.pePrw .53535+Cf 
END IF 
IF Ihi=1]. THEN 
Sra=.88-( .24/(4.+Prwf )) 
SP DS5 935345 6S eo ent rete 
Ome=(5.+.015*Reif *Sra*Prwf*Srb ) 
END he 
Dehn 2 eRe ns 
Epsi=( 1282 *LGl( Re: j)-1.64) (2) 
Pok |=) 343 34*eps i 
Pak ei. feces au =) 
Polt=(Epsi/8)*Rei*Prw 
ser Maisie <2 3Ulsp cies) psa riem lia = | i) 
Ome=Fol/Fac2 


104 


ene .F 

Hi=rw/D1*L1* ome 

hee = -O THEN bitd 249) 
eee le ot) 

Per tO, Fie Sie-Od + .5 
MlHtCHierl. Fea ed |» 
mee te (Oi +s * 

feet -)) )+P (by +e 
Meyer Ve CuUtAL sy) °.S 
meee Nianht Mie o7¢Ml «Ll ) 
mee--Niann(Metl ls (Molec } 


wn 


Ot=O/(P1l*Oiet(bLe+Li«Felt+tLieFer '#H1 ) 


Ir Thi=8 THEN 


Sie— How/PNMuut Tavgo*Ot:?”. 


If mbsi (Che -Cho/Cfe. .0Q01 


ENO 
ENO oLr 
few thi=! THEN 

peurcec=eNrrwt lavgtDt } 


Reirfc=Vwe01*FNRhow( TavgtDt }/FNMuw( Tavgt+Dt ) 
eee ot (Prure—-eprcut )/Prwute)>.00} 


Pru -Seruretecwt 7. 
Rerf=(Reife+Reif )/2. 
GCTU 2451 
END LF 
BNO. iF 
ie lfto 


meee = 120-111 )7LO6(%i Tsteam-T11)/( Tsteam-T2o)) 


Uo=Q/iLmtd*Pl*DorL ) 


pee | 7lo-Yo*L/(Di*tLtlLi*Fel +lLo*#Fec )*H: »-Rm) 


Tef=Qp/Ho 

Cosc=FNCpwl! Tconrlsteam)+#.§S) 
Hf g=FNHfgi(Tsteam) 
Two=Tsteam-9p/Ho 
Tfhilm=Tsteam/3+Twor2/3 
Kf=FNKw( Thilm) 

Rhof=FNRhow( Tfilm) 

Mu f=FNMuw( Tfilm) 


ae. Gol *if*(Rhof 2*9.81*Hfa/(Muf*00#0p))*.3333 
Bim@iss=-7eo*tht 3*9.98)]*Hfg*Rhof 2c/(Muf*Do*lcf))*.cS 


Alpl=.728*#Ho/Hnuss 
fem d-l=lfilm 
Qpa(J~-1)=Qp 
Tompetyo 49000 


| 4 
THEN 


O'S 


OR ABS( (Rerte=Reit )/Reife), .vgl 


“557 «=p 
2560 S.=S<4k 
2563 Sy=Syev 
2665 645-5950 2 
[SBS Saysdnyteey 





set. = Ul-Sug 


cS75 OQloss=017( 100-25) [stean-(ro0om) 9s) 1Gkee Uri a re. 
S76 Htc=FNHf* feon 

~584 Mdv=v 

“S87 Bp=(Bvol+!10@:;°2/5.76 


25390) HSc=lase*t ict san-deerw 

e595 Mdve=((bp-U less) -Mdv* tise, Ai 

2-596 IF ABS((Mevatdve Maeve a ae | 
2599 Mdv=(Mdv+Mdvc j#.5 | 
e602 GOTO 2533 | 
SoS ee aa 
~608 Mdv='MdveMdve 25 


-611t VUg=FNUvst'Tsteam) 
~614 Uv=emMdv#¥g/Ayz 
e620 =. So2\ Dour *Ato 7/0 UN 2k eC Steam iwome 


2625" NU=Horioe, kt 

e626  ReteVyslmhot *De/Mut 

e629 Nr=Nu/Ret*.5 

550 Hfuj3=.96*(9.8t*Hro/ict )°22ekKt .88Uy” Seiko 25 Go oe eee 
eboS LF lov=l een 


2645 A ONT 

fe ae Pein USTiGS S72 Vu Rei Hi Uo Hfujy(OT) 
Hnui Qo" "" 

-652 PRINT USING “SX ,2.00,1X ,3(MZ.30E ,1X) 3k ,2(MZ.30E ,5X)° 3 UW Rea 

sinleie 

Soe: PRINT 

Pass Ee UN US 1 NG aS Krs my, Ho Q Tet NuRe F 
Slate Cee 

21) Ve PRINT USING “SX .2.00,1X ,c(MZ.30E ,1X),2X ,3D.00 2X ,3(MZ. 30E | 

ee Nie e ass 

file s\e PRINT 


oocg eee 

26593 [Ff lev=) THEN 

-662 eo Ugale “ealent 

266! PRINT USING “11X ,00,2xX ,Z.00,1X ,3(MO.30E ,1X) ,2(30.00,1x% ) O 7000 eee 
He On .icf. [stmam send 

26p iL oe 

2668 PRINT USING “11X ,00,4X ,Z.00 ,2X ,2(M0.30E ,2X) ,Z.30E ,3X ,50.00 ,2X% -SU eee 
We Ho Ope tet sis team 

palsy qe) ale 

-674 ENO IF 
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-j UI) @) OF 
& toe =) 
Si) ges fey yy 


3] 


tc 


ius Op 


4 


es Vee Vis) re Feral el oy Tp eg es 
J =d =. J] 


mint! mm — 


J 
ae 


‘ —- —— —_ — ed 


r 


Oo a) Cis 


oui] 25 ff fh ie) [ates 


Go GW mM as wm 


OM — G7) Ul & C1 Cl me fat: 


FIORI hehe fae eri ee a a 3 


ty) WN WH UW Of Cd & 


wo G& « aw 


‘ Poe aes so=1 be iw seine eUTPUT SFicep;Cp Ho 


@ m=) THEN 
Ce ee ee OD) Per a ee N OUTPUT “er leps a ,Ho 
Mee ee CaANcCe “(COL RISE? l=¥, c=N" Ltr 
See te =| HEN Soto Coed 
Siam 
Mme tees ORE THIS CATA SET ly, 8=N)?" Oks 


Pee; plese vol. bamp etaven.)) .12.Pvaplarvaod emt Ue) 
fieeerlo)lte lon 


Bese 
[=o -1 
on aes 
Beer 
PHPUT WILL THESE SE RNOTHER RUN = 1=4 9=N)?" Go_on 
Nrun= J 
PP 50 orn @ (HEN Fepeat 
Boe 
Peoria $b helsRepeat 
eNO iF 
Beeoe=! THEN 220s 
IF ITwrl=a THEN 
pect oN Berle 7D « 
1job=i 
iwd= | 
Lm=2 
VESIGN®OPi le TO DO files 
SOT 1156 
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